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Rapid  progress  in  the  synthesis  and  processing  of  materials  with  structure  on  nanometer  length 
scales  has  created  a  demand  for  greater  scientific  understanding  of  thermal  transport  in  nanoscale 
devices,  individual  nanostructures,  and  nanostructured  materials.  This  review  emphasizes 
developments  in  experiment,  theory,  and  computation  that  have  occurred  in  the  past  ten  years  and 
summarizes  the  present  status  of  the  field.  Interfaces  between  materials  become  increasingly 
important  on  small  length  scales.  The  thermal  conductance  of  many  solid- solid  interfaces  have  been 
studied  experimentally  but  the  range  of  observed  interface  properties  is  much  smaller  than  predicted 
by  simple  theory.  Classical  molecular  dynamics  simulations  are  emerging  as  a  powerful  tool  for 
calculations  of  thermal  conductance  and  phonon  scattering,  and  may  provide  for  a  lively  interplay 
of  experiment  and  theory  in  the  near  term.  Fundamental  issues  remain  concerning  the  correct 
definitions  of  temperature  in  nonequilibrium  nanoscale  systems.  Modern  Si  microelectronics  are 
now  firmly  in  the  nanoscale  regime — experiments  have  demonstrated  that  the  close  proximity  of 
interfaces  and  the  extremely  small  volume  of  heat  dissipation  strongly  modifies  thermal  transport, 
thereby  aggravating  problems  of  thermal  management.  Microelectronic  devices  are  too  large  to 
yield  to  atomic-level  simulation  in  the  foreseeable  future  and,  therefore,  calculations  of  thermal 
transport  must  rely  on  solutions  of  the  Boltzmann  transport  equation;  microscopic  phonon  scattering 
rates  needed  for  predictive  models  are,  even  for  Si,  poorly  known.  Low-dimensional  nanostructures, 
such  as  carbon  nanotubes,  are  predicted  to  have  novel  transport  properties;  the  first  quantitative 
experiments  of  the  thermal  conductivity  of  nanotubes  have  recently  been  achieved  using 
microfabricated  measurement  systems.  Nanoscale  porosity  decreases  the  permittivity  of  amorphous 
dielectrics  but  porosity  also  strongly  decreases  the  thermal  conductivity.  The  promise  of  improved 
thermoelectric  materials  and  problems  of  thermal  management  of  optoelectronic  devices  have 
stimulated  extensive  studies  of  semiconductor  superlattices;  agreement  between  experiment  and 
theory  is  generally  poor.  Advances  in  measurement  methods,  e.g.,  the  3co  method,  time-domain 
thermoreflectance,  sources  of  coherent  phonons,  microfabricated  test  structures,  and  the  scanning 
thermal  microscope,  are  enabling  new  capabilities  for  nanoscale  thermal  metrology.  ©  2003 
American  Institute  of  Physics.  [DOI:  10.1063/1.1524305] 
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L  INTRODUCTION 

Recent  advances  in  synthesis,  processing,  and  mi¬ 
croanalysis  are  enabling  the  routine  production  of  well- 
characterized  materials  with  structure  that  varies  on  the 
length  scale  of  several  nanometers.  Examples  are  semicon¬ 
ductor  quantum  dots  and  superlattices,  polymer  nanocom¬ 
posites,  multilayer  coatings,  microelectronic  and  optoelec¬ 
tronic  devices,  and  microelectromechanical  sensors.  Many  of 
these  nanoscale  structures  already  have  important  commer¬ 
cial  applications,  while  others  are  studied  scientifically.  The 
most  important  devices,  in  terms  of  worldwide  sales,  are  the 
integrated  circuits  used  in  computer  circuits  and  memories. 
Device  engineers  have  been  astonishingly  clever  at  maintain¬ 
ing  Moore’s  Law,  which  states  that  the  individual  device 
components  get  smaller  every  year.  Soon  the  field-effect- 
transistors  (FET)  will  have  a  channel  length  of  100  nm,  and 
50  nm  devices  are  projected  in  the  near  future.  In  most  cur¬ 
rent  and  envisioned  applications  of  nanostructures,  thermal 
management  is  a  serious  issue.  In  some  devices,  such  as 
computer  processors  or  semiconductor  lasers,  one  wants  to 
get  the  heat  away  as  efficiently  as  possible — these  systems 
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need  high  thermal  conductivity.  In  others,  such  as  thermal 
barriers  or  thermoelectric  materials  used  for  solid-state  re¬ 
frigeration,  one  wants  the  thermal  conductivity  to  be  as  small 
as  possible. 

Our  discussion  emphasizes  thermal  transport  in  nonme- 
tallic  systems,  in  which  heat  is  transported  by  phonons. 
Phonons  have  a  wide  variation  in  frequency,  and  an  even 
larger  variation  in  their  mean-free-paths  (mfps).  However, 
the  bulk  of  the  heat  is  often  carried  by  phonons  of  large  wave 
vector,  and  they  have  mfps  at  room  temperature  of  1-100 
nm.  So  in  many  systems  of  current  interest,  the  scale  of  the 
microstructure  is  the  same  scale  as  the  mfp  of  the  phonons, 
and  sometimes  comparable  to  the  phonon  wavelength.  This 
clearly  necessitates  an  understanding  of  heat  transport  be¬ 
yond  that  achievable  at  the  continuum  level.  So  far,  no  ana¬ 
lytical  theories  have  adequately  treated  the  wave  nature  of 
phonons  that  will  be  manifest  at  these  length  scales.  More¬ 
over,  this  equality  in  length  scales  raises  important  concep¬ 
tual  questions,  among  which  is  how  to  define  temperature  on 
the  length  scale  of  the  phonon  wavelength  and  phonon  mfp. 

This  review  will  summarize  the  present  understanding  of 
thermal  transport  in  nanostructures.  The  experimental  diffi¬ 
culties  begin  with  the  severely  limited  capabilities  for  mea¬ 
suring  thermal  transport  in  increasingly  small  systems.  A  re¬ 
cent  review  summarizes  the  experimental  techniques  that  are 
being  used  to  probe  thermal  conduction  at  a  submicron 
level.1  These  methods  are  also  discussed  here.  Currently,  the 
only  method  with  nanometer- scale  spatial  resolution  that 
may  be  applied  to  nanometer  structures  is  based  on  the 
atomic  force  microscope  (AFM).  However,  AFM  measure¬ 
ments  are  only  beginning  to  provide  an  understanding  of 
thermal  transport  in  nanoscale  structures;  the  experimental 
field  is  certainly  in  its  infancy. 

The  theory  and  simulation  of  nanoscale  thermal  trans¬ 
port  is  at  a  similarly  immature  stage.  The  recognition  that  a 
more  detailed  theoretical  understanding  of  heat  transport 
than  achievable  using  Fourier’s  Law  of  heat  conduction,  Q 
=  —  ac V  T  where  Q  is  the  heat  flux,  k  the  thermal  conductiv¬ 
ity  and  VT  the  temperature  gradient,  has  lead  to  the  devel¬ 
opment  of  approaches  based  on  the  numerical  solutions  of 
the  Boltzmann  transport  equation  (BTE)  and  to  atomic-level 
simulations  of  thermal  conductivity.  Each  still  requires  fur¬ 
ther  methodological  developments  and  the  first  systematic 
applications  of  these  methods  to  heat  transport  issues  are  still 
being  refined. 

Since  interfaces  play  a  critical  role  in  nanoscale  thermal 
transport,  we  begin  our  discussion  with  a  review  of  experi¬ 
ments  and  theory  on  the  transport  of  heat  between  two  dif¬ 
ferent  materials  or  between  twin  or  twist  boundaries  of  the 
same  material.  In  Sec.  Ill,  we  focus  on  thermal  transport  in 
microfabricated  nanostructures  and  synthetic  nanostructures 
of  reduced  dimension.  The  related  topic  of  heat  transport  in 
nanostructured  thin  film  materials,  including  superlattices,  is 
addressed  in  Sec.  IV.  We  explore  issues  associated  with  me¬ 
trology  and  measurement  science  in  Sec.  V. 
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II.  HEAT  TRANSPORT  AND  PHONON  DYNAMICS  AT 
INDIVIDUAL  INTERFACES 

The  idea  that  an  interface  should  produce  a  thermal  re¬ 
sistance  is  intuitively  appealing,  since  an  interface  consti¬ 
tutes  an  interruption  in  the  regular  crystalline  lattice  on 
which  phonons  propagate.  For  an  interface  between  dissimi¬ 
lar  materials,  the  different  densities  and  sound  speeds  result 
in  a  mismatch  in  the  acoustic  impedances;  this  is  directly 
analogous  to  the  mismatch  in  the  refractive  indices  of  two 
optically  different  materials.  The  effect  that  this  impedance 
mismatch  has  on  phonon  transmission  is  captured  in  the 
acoustic-mismatch  (AM)  model.  By  assuming  that  no  scat¬ 
tering  takes  place  at  the  interface,  and  by  imposing  appropri¬ 
ate  stress  and  displacement  boundary  conditions  at  the  inter¬ 
face,  the  AM  model  gives  the  transmission  coefficient  tAB  for 
phonon  energy  in  material  A  incident  normal  to  the  interface 
with  material  B  as: 


^ZAZB 

tAB~(ZA  +  ZB)2’ 


(1) 


where  Z=  pc  is  the  acoustic  impedance,  with  c  and  p  being 
the  speed  of  sound  and  mass  density,  respectively.2,3  In  the 
AM  model,  the  interface  has  no  intrinsic  properties  but 
merely  joins  the  two  grains.  As  a  consequence  a  twist  grain 
boundary,  in  which  the  normals  on  the  two  sides  of  the  GB 
are  crystallographically  identical,  would  incorrectly  be  pre¬ 
dicted  to  have  a  transmission  coefficient  of  one. 

By  contrast,  in  the  diffuse  mismatch  (DM)  model  it  is 
assumed  that  all  phonons  striking  an  interface  loses  memory 
of  where  they  came  from.  Then,  the  probability  of  being 
scattered  to  one  side  of  the  interface  or  the  other  is  simply 
proportional  to  the  phonon  density  of  states.  Thus,  as  in  the 
AM  model,  the  fraction  of  energy  transmitted  is  independent 
of  the  structure  of  the  interface  itself.3 

While  the  above  two  models  provide  useful  reference 
calculations  against  which  to  compare  experimental  results, 
and  for  many  cases  give  rather  similar  predictions,  neither 
captures  the  complexity  of  the  interaction  between  phonons 
and  real  interfaces. 

As  we  will  see  in  this  section,  the  interactions  of 
phonons  with  a  single  interface  still  offers  significant  chal¬ 
lenges  to  both  experiments  and  theory/simulation.  We  first 
review  experiments  on  thermal  transport  across  individual 
interfaces,  including  both  the  classic  Kapitza  problem  of  a 
superfluid-helium/solid  interface  and  model  solid/solid  inter¬ 
faces.  We  point  out  that  spectral  methods,  involving  phonons 
of  well-defined  frequencies  and  wave  vectors,  offer  the 
promise  of  providing  insights  into  thermal  transport  not  ac¬ 
cessible  to  more  traditional  experiments.  We  summarize  re¬ 
cent  atomic-level  simulations  of  thermal  conductivity 
through  individual  interfaces  and  identify  a  need  for  simula¬ 
tions  of  the  interaction  of  individual  phonons  with  interfaces. 
Finally,  we  point  out  that  the  simulations  of  thermal  conduc¬ 
tivity  raise  fundamental  thermodynamics  issues  concerning 
the  definition  of  temperature  of  nonequilibrium  nanoscale 
systems. 


A.  Thermal  conductance  of  model  interfaces 

We  consider  the  thermal  conductance  of  an  isolated  in¬ 
terface,  i.e.,  an  interface  that  is  separated  from  other  inter¬ 
faces  by  a  distance  that  is  large  compared  to  the  mean-free- 
path  of  the  lattice  vibrations  that  dominate  heat  transport  in 
the  material.  In  this  limit,  we  can  ignore  coherent  superposi¬ 
tion  of  lattice  waves  reflected  or  transmitted  by  adjacent  in¬ 
terfaces.  Swartz  and  Pohl3  give  an  exhaustive  review  of  ther¬ 
mal  boundary  resistance,  the  inverse  of  interface  thermal 
conductance,  through  1988. 

The  existence  of  a  thermal  boundary  resistance  between 
a  solid  and  superfluid  helium  was  first  detected  by  Kapitza  in 
Ref.  4.  When  the  discovery  was  first  made,  it  was  thought 
that  this  resistance  was  due  to  some  special  property  of  su¬ 
perfluid  helium.  Only  later  was  it  realized  that  such  a  Kaptiza 
resistance,  RK ,  exists  at  the  interface  between  any  pair  of 
dissimilar  materials.  Khalatnikov2  developed  the  acoustic- 
mismatch  model  to  explain  this  thermal  resistance.  The 
acoustic  mismatch  is  extremely  large  at  an  interface  between 
helium  and  most  ordinary  solids  because  the  density  and 
sound  velocity  of  helium  are  very  much  less  than  for  other 
materials.  While  Khalatnikov’ s  theory  was  able  to  explain 
the  existence  of  the  Kapitza  resistance,  the  value  of  RK  that 
was  calculated  was  between  10  and  100  times  larger  than 
typical  experimental  values  measured  at  1-2  K.  This  indi¬ 
cates  that  phonons  are  able  to  cross  the  interface  substan¬ 
tially  more  easily  than  expected. 

Remarkably,  the  reason  for  this  so-called  “anomalous 
transmission”  at  these  low  temperatures  has  still  not  been 
established.  It  was  proposed  by  Challis,  Dransfeld,  and 
Wilks5  that  the  transmission  was  increased  due  to  the  exis¬ 
tence  of  a  compressed  layer  of  helium  on  the  surface  of  the 
solid.  This  layer  is  formed  because  of  the  van  der  Waals 
attraction  of  helium  atoms  to  a  solid  surface.  An  important 
experiment  was  performed  by  Weber  et  al6  They  measured 
the  transmission  of  phonons  across  the  interface  between  an 
alkali-halide  crystal  and  helium.  Very  clean  surfaces  were 
produced  by  cleaving  the  alkali-halide  crystals  at  low  tem¬ 
perature.  For  these  surfaces  no  anomalous  transmission  was 
observed,  indicating  that  the  excess  transmission  must  be 
associated  with  some  sort  of  surface  contamination  or  de¬ 
fects.  Olson  and  Pohl7  verified  that  the  anomalous  transmis¬ 
sion  between  Si  and  liquid  helium  directly  correlates  with  an 
increase  in  the  diffuse  phonon  scattering. 

Since  the  work  of  Weber  et  al.  6  a  number  of  attempts 
have  been  made  to  explain  how  surface  contamination  can 
have  such  a  large  effect  on  the  transmission.  Much  of  this 
work  is  reviewed  in  Refs.  3,  8,  and  9.  Despite  these  efforts, 
no  generally  accepted  theory  has  emerged.  The  resolution  of 
this  problem  requires  more  detailed  experimental  work  under 
controlled  conditions  with  well-characterized  surfaces  pre¬ 
pared  and  studied  under  ultrahigh  vacuum. 

Of  course,  while  not  directly  relevant  to  technical  prob¬ 
lems  of  thermal  management  in  nanostructures,  the  work  on 
superfluid  helium  does  leave  open  the  possibility  that  the 
transmission  of  phonons  across  the  interface  between  two 
solids  may  be  substantially  affected  by  a  small  number  of 
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FIG.  1.  Comparison  of  selected  data  for  interface  thermal  conductance:  (i) 
individual  interfaces  measured  by  picosecond  thermoreflectance  (Ref.  10), 
A1/A1203  (open  diamonds)  and  Pb/c-C  (filled  triangles);  (ii)  conductance  of 
the  «-GeSbTe2.5  /ZnS  interface  (open  triangles)  from  a  multilayer  sample 
(Ref.  14),  (iii)  series  conductance  of  the  top  and  bottom  interfaces  of  metal  - 
Si02-silicon  structures  (filled  circles  labeled  by  MOS)  (Refs.  12  and  13). 
The  solid  line  is  the  calculated  diffuse  mismatch  conductance  of  A1/A1203 
using  the  Debye  model;  the  dashed  line  is  the  theoretical  prediction  for 
A1/A1203  using  a  lattice-dynamical  calculation  of  a  model  fee  interface 
(Ref.  10).  The  right  axis  gives  the  thickness  of  a  film  with  A 
=  1  W  m~ 1  K“ 1  that  has  the  thermal  conductance  corresponding  to  the  left 
axis. 


surface  defects.  The  extent  to  which  this  occurs  has  not  yet 
been  established. 

Some  recent  data  for  the  thermal  conductances  of  model 
solid/solid  interfaces  are  shown  in  Fig.  1.  Surprisingly,  near 
room  temperature,  the  highest  thermal  conductance  for 
metal/dielectric  interfaces  measured  by  picosecond 
thermoreflectance,10  Ti/Al203  and  A1/A1203 ,  is  only  a  factor 
of  ~5  larger  than  the  lowest  conductance,  Pb/diamond. 
Other  data  for  metal/dielectric  interfaces,3,11-13  for  the  most 
part,  fall  between  these  two  extremes. 

B.  Heat  pulse  and  coherent  phonon  experiments 

Phonon  transmission  across  a  single  interface  has  been 
studied  using  both  coherent  and  incoherent  phonons.  To 
study  incoherent  phonons  a  pulse  of  energy  is  deposited  into 
the  material  on  one  side  of  the  interface.  This  generates  in¬ 
coherent  phonons  with  a  frequency  spectrum  that  is  deter¬ 
mined  by  the  pulse  energy  and  the  ambient  temperature. 
These  phonons  propagate  across  the  interface  and  are  then 
detected  in  the  material  on  the  other  side.  This  type  of  ex¬ 
periment  has  been  performed  for  interfaces  between  solids 
and  liquid  helium  (see,  for  example,  Wichert  et  al. 15)  and  at 
solid-solid  interfaces  (for  a  review,  see  Ref.  16).  Usually,  the 


FIG.  2.  (Color)  Top:  experimentally  observed  phonon  focusing  image  in 
twist-bonded  Si  with  a  twist  angle  of  40  degrees  (Ref.  17).  Each  of  the  two 
bonded  crystals  are  3  mm  thick,  the  width  of  the  image  corresponds  to  2 
mm.  Bottom:  simulated  image.  Arrows  indicate  the  crystal  orientation  rela¬ 
tive  to  the  phonon  source.  Phonon  intensities  are  normalized  to  the  maxi¬ 
mum  value  in  each  image. 


experiments  are  performed  at  low  temperatures  so  that  the 
phonon  mean  free  path  is  long;  i.e.,  the  phonons  travel  bal- 
listically  in  the  two  media.  As  a  result,  it  is  possible  to  study 
separately  the  transmission  of  phonons  of  different  polariza¬ 
tion  and  to  observe  the  refraction  that  occurs  when  the 
phonons  cross  the  interface.  This  refraction,  combined  with 
the  effects  of  phonon  focusing  lead  to  an  interesting  angular 
variation  of  the  phonon  flux  emerging  from  the  interface.16 
As  an  example,  in  Fig.  2  we  show  the  phonon  flux  emerging 
from  the  interface  between  two  twist-bonded  Si  crystals.17 
The  faces  of  these  crystals  both  have  a  (100)  orientation,  but 
one  crystal  is  rotated  by  40°  around  the  [100]  axis  relative  to 
the  other. 

To  study  the  propagation  of  coherent  phonons  across  an 
interface,  the  picosecond  ultrasonic  technique  can  be  used.  A 
short  light  pulse  is  absorbed  in  a  material,  and  because  of 
thermal  expansion  a  short  strain  pulse  (coherent  phonon)  is 
generated.  This  strain  pulse  is  detected  after  it  has  been  re¬ 
flected  or  transmitted  across  the  interface  to  another  material. 
In  this  way  the  transmission  of  phonons  of  frequency  up  to 
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about  1  THz  can  be  studied.  In  one  experiment  of  this  type, 
coherent  phonon  transmission  at  a  bonded  interface  between 
Si02  and  Si  was  studied.18,19  It  was  possible  to  use  the 
acoustic  transmission  as  a  measure  of  the  quality  of  the  bond, 
and  to  see  how  the  bonding  improved  with  thermal  anneal¬ 
ing.  In  another  experiment,20  a  graded  interface  between 
films  of  gold  and  aluminum  was  studied.  The  interface  width 
£  was  ~  400  A .  For  a  diffuse  interface  the  phonon  transmis¬ 
sion  is  expected  to  depend  on  the  value  of  k£,  where  k  is  the 
phonon  wave  number.  For  1,  the  transmission  should  be 
larger  because  the  acoustic  properties  vary  smoothly  on  the 
scale  of  the  phonon  wavelength.  The  experiment  demon¬ 
strated  this  effect. 

Quantitative  measurements  with  coherent  phonons  have 
many  potential  applications  for  the  characterization  of  inter¬ 
faces.  In  picosecond  ultrasonics,  the  light  creates  sound 
pulses,  i.e.,  a  continuum  of  modes  for  which  the  relevant 
length  scale  is  defined  by  the  penetration  depth  of  the  light 
and  other  parameters  of  the  overlayer.  Periodic  layered  me¬ 
dia  offer  an  alternative  method  for  generating  coherent 
phonons.  Here,  phase-matching  considerations  show  that 
light  can  couple  to  modes  of  wave  vectors  q^lirn/d  where 
n  is  an  integer  and  d  is  the  period  of  the  structure.  Using 
femtosecond  laser  sources,  this  approach  has  been  exten¬ 
sively  applied  to  GaAs-Al^Gai  _xAs  (Ref.  21)  and 
In^Ga^^N-GaN  (Refs.  22  and  23)  superlattices  to  generate 
coherent  acoustic  modes  of  frequencies  up  to  ~  1  THz. 

C.  Molecular-dynamics  simulation  of  heat  transport 
and  phonon  dynamics 

Computational  approaches  to  heat-transfer  problems 
span  the  range  from  numerical  solutions  of  Fourier’s  law24  to 
calculations  based  on  the  Boltzmann  transport  equation 
(BTE)25,26  to  atomic-level  simulations.  In  the  Fourier  and 
BTE  approaches,  the  physics  of  heat  transfer  and  phonon 
scattering  are  incorporated  into  the  calculations  in  an  explicit 
manner;  thus  for  a  reliable  calculation,  a  fairly  sophisticated 
understanding  of  the  fundamental  phonon  processes  is  re¬ 
quired.  While  in  some  cases  (e.g.,  anharmonic  phonon  ef¬ 
fects,  isotopic  defects,  and  point  defects)  understanding  is 
well-developed,  in  others  (e.g.,  interfaces)  current  under¬ 
standing  is  poor.  Molecular-dynamics  (MD)  simulation,  by 
contrast,  merely  involves  the  integration  in  time  of  Newton’s 
equations  of  motion  for  an  ensemble  of  atoms  interacting 
with  each  other  through  a,  usually  empirical,  interatomic 
potential.27,28  Because  the  formalism  of  the  MD  approach 
does  not  require  any  a  priori  understanding  of  heat  transport, 
it  is  ideal  for  investigating  the  fundamental  heat-transfer 
mechanisms  themselves.  However,  MD  does  have  the  sig¬ 
nificant  limitation  of  being  entirely  classical,  with  each  vi¬ 
brational  mode  equally  excited;  thus  it  is  only  rigorously 
applicable  to  solids  above  the  Debye  temperature.  Moreover, 
since  electrons  are  not  included  in  an  atomistic  model,  it  is 
not  possible  to  simulate  electrical  conductors  or  the 
electron-phonon  interactions  present  in  many  semiconduc¬ 
tors.  In  this  subsection  we  briefly  review  recent  atomic-level 
simulations  of  thermal  transport  at  interfaces,  and  point  out 
the  opportunities  for  gaining  important  new  insights  from  the 
simulation  of  individual  phonon  processes. 
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FIG.  3.  Thermal  profile  calculated  by  MD  for  a  X  =  13  twin  boundary  in  Si, 
after  Ref.  31.  The  sample  is  25  nm  along  the  x  axis  and  contains  840  atoms. 
Total  simulation  time  is  1  ns. 


The  two  most  commonly  used  approaches  to  the  MD 
simulation  of  the  thermal  conductivity  are  the  Green- Kubo 
approach  in  which  the  equilibrium  fluctuations  in  the  heat 
current  are  analyzed,  and  the  “direct  method,”  which  mimics 
experiment  by  imposing  a  temperature  gradient  on  the  sys¬ 
tem  and  determining  the  thermal  conductivity  from  Fourier’s 
Law  (for  a  recent  review  of  these  methods  see  Ref.  29). 
While  simple  in  principle,  in  practice  each  has  significant 
complications  associated  with  system-size  effects  and  the 
relatively  long  simulation  times  required.  Moreover,  the  di¬ 
rect  method  has  the  additional  complication  of  a  potential 
nonlinear  response  of  the  system  associated  with  the  very 
large  temperature  gradients  imposed  (typically  in  excess  of 
109  Km-1).  A  recent  systematic  comparison  has  shown  that, 
when  extrapolated  to  infinite  system  size,  the  Green-Kubo 
and  direct  methods  give  consistent  values  for  the  thermal 
conductivity,  albeit  each  with  error  bars  of  20%  or  more.29 
Despite  the  additional  complications,  it  was  concluded  that 
the  direct  method  is  usually  preferable,  particularly  when 
applied  to  interfacial  systems,  as  it  has  the  significant  advan¬ 
tage  of  providing  a  direct  measure  of  the  interfacial  resis¬ 
tance,  Rk  ,  from  the  temperature  drop  across  the  interface. 
(Unfortunately,  however,  as  discussed  in  the  next  subsection, 
the  conceptual  issues  associated  with  defining  the  tempera¬ 
ture  in  an  inhomogeneous  system  have  yet  to  be  fully  re¬ 
solved.) 

There  have  been  relatively  few  atomic-level  simulations 
of  heat  transport  at  solid- solid  interfaces.  Pickett  et  al .30  ex¬ 
plored  the  transmission  of  phonon  energy  at  an  interface  be¬ 
tween  two  diamond-structured  materials  differing  only  in 
their  masses.  Maiti  et  al .31  used  the  direct  method  to  perform 
the  first  simulations  of  thermal  transport  through  grain 
boundaries.  Figure  3  shows  the  temperature  profile  through 
the  simulation  cell.  The  thermal  conductance,  GK,  of  the 
interface  can  be  calculated  from  the  temperature  drop  A  T, 
determined  from  the  figure,  and  from  the  heat-current  density 
Q  passing  through  the  system,  determined  from  the  simula¬ 
tion:  GK=  Q/AT.  For  two  different  symmetric  tilt  grain 
boundaries  in  Si,  they  obtained  G^=0.8  and  GK 
=  0.9  GW  m-2  K-1. 
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Recently,  Schelling  et  al.32  determined  the  thermal  con¬ 
ductance  for  the  (111)  231  twist  GB,  which  has  a  rather  low 
energy,  and  for  the  (100)  229  twist  GB,  which  has  a  rather 
high  energy,  for  the  same  interatomic  description  of  Si  and 
obtained  values  of  1.53  and  0.85  GW m-2  K-1,  respectively. 
In  a  number  of  studies  of  the  GBs  in  Si,  Keblinski  et  al  33 
showed  that  the  (100)  229  GB  is,  from  a  structural  point  of 
view,  representative  of  the  many  high-energy  GBs  in  Si.  To¬ 
gether  with  the  very  similar  values  for  the  two  tilt  GBs,  this 
suggests  that  G^0.8-0.9GWm“2K_1  may  be  taken  as 
representative  of  most  GBs  in  Si.  The  much  higher  value  for 
conductance  of  the  (111)  231  is  reasonable  given  that  the 
GBs  on  (111)  plane  have  the  lowest  energy  of  all  the  high- 
angle  GBs  in  Si.33  In  comparison  with  the  experimental  re¬ 
sults  show  in  Fig.  1,  these  values  for  the  conductance  are 
rather  high.  One  possible  reason  for  this  is  that  all  of  the 
experimental  systems  are  interfaces  between  dissimilar  ma¬ 
terials  which,  due  to  the  acoustic-impedance  mismatch,  can 
be  expected  to  have  a  lower  conductance  than  a  grain  bound¬ 
ary.  A  second  possible  reason  could  be  that  the  simulated 
grain  boundaries  are  defect- free  and  atomically  flat,  thus  pre¬ 
sumably  minimizing  the  thermal  resistance. 

It  is  useful  to  analyze  the  thermal  resistance  of  an  inter¬ 
face  in  terms  of  the  length  of  perfect  crystal  that  would  pro¬ 
vide  an  equivalent  thermal  resistance.  This  led  Nan  et  al34  to 
define  the  Kapitza  length,  Ik=k/Gk  ;  typical  values  for  Si 
(ac—  100  W  m-1  K-1  and  GK~  lGWm"2K_1)  give  lK 
—  100  nm  as  the  thickness  of  perfect  crystal  with  the  same 
thermal  resistance  as  the  interface. 

Figure  3  shows  regions  of  nonlinear  temperature  change 
far  from  the  GB.  These  nonlinear  regions  are  close  to  the 
heat  source  and  heat  sink  in  the  system  and  arise  from  the 
interface  between  the  regions  that  are  heated  and  cooled  and 
the  rest  of  the  system.  The  interpretation  of  these  nonlinear 
regions  can  draw  from  the  results  of  phonon  transport 
theory,25,26  which  has  connections  to  the  kinetic  theory  of 
molecular  transport35,36  and  the  equation  of  photon  radiative 
transfer.37,38  These  increased  temperature  gradients  resemble 
those  predicted  by  the  Boltzmann  transport  equations  for 
gases,  photons,  and  phonons. 

Isothermal  boundary  conditions  enforce  perfect  emission 
and  absorption  of  phonons  from  the  boundary,  i.e.,  a  phonon 
blackbody,  and  ensure  equilibrium  conditions  for  the  emitted 
phonons.  This  has  an  influence  comparable  to  an  ensured 
phonon  scattering  event  at  the  boundary.  Phonons  traveling 
away  from  the  boundary  have  a  distribution  characteristic  of 
the  boundary  temperature,  which  differs  from  the  distribution 
of  phonons  traveling  in  the  opposite  direction  more  than  the 
difference  found  deep  within  the  medium.  This  results  in  an 
increased  departure  from  equilibrium  locally  within  the  pho¬ 
non  system  and  an  increased  gradient  in  the  directionally 
averaged  phonon  energy.  This  is  called  the  temperature  slip 
phenomenon,  and  is  well  documented  theoretically  for  gas 
molecules,  photons,  and  phonons.  It  occurs  within  a  length 
scale  of  the  boundary  comparable  with  the  mean  free  path. 
The  magnitude  of  the  temperature  slip,  specifically  the  dif¬ 
ference  between  the  boundary  temperature  and  the  extrapo¬ 
lation  of  the  linear  temperature  distribution  deep  within  the 
medium,  is  approximately  the  product  of  the  heat  flux  and 
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FIG.  4.  Interaction  of  a  wavepacket  of  longitudinal  acoustic  phonons  [o> 
^8.7  THz,  kz^ 0.62(2 ir/a)]  with  an  (001)  epitaxial  interface  between  two 
diamond-structured  lattices  differing  only  in  the  mass  of  the  atoms  (mass 
=  1  on  the  right  of  the  interface,  mass  =  4  on  the  left  of  the  interface);  37% 
of  the  energy  in  the  incoming  acoustic  phonon  package  is  transmitted  in  the 
form  of  a  slow-moving  optical  phonon,  while  the  remaining  63%  is  reflected 
with  no  change  in  frequency  or  wavelength  (Ref.  42). 


the  ratio  of  the  mean  free  path  and  the  bulk  thermal  conduc¬ 
tivity.  This  argument  is  based  on  the  Rosseland  diffusion 
approximation.39  The  near-boundary  predictions  of  atomistic 
simulations  should  be  compared  with  the  solutions  to  the 
phonon  transport  equation  in  the  vicinity  of  boundaries,  or  at 
least  to  the  closed-form  expressions  for  the  temperature  slip 
phenomenon. 

Although  the  value  for  the  thermal  conductance  of  a  par¬ 
ticular  interface  in  a  particular  material  may  be  a  useful 
quantity  to  know,  it  provides  no  mechanistic  information. 
The  explicit  relationship  between  the  temperature-dependent 
thermal  conductance,  GK(T ),  and  the  phonon  transmission 
coefficient  at  an  interface  has  been  given  by  Young  and 
Maris:40 

Gk(T)=  f  da>C(co,T)(vz(co))(t(co))  (2) 

where  C(a),T)  is  the  specific  heat,  which  can  be  determined 
from  the  phonon  density  of  states  and  the  Bose-Einstein 
distribution  function,  vz((o)  is  the  component  of  the  phase 
velocity  of  the  phonon  normal  to  the  interface  and  t(oj)  is 
the  transmission  coefficient.  The  (  )  denotes  an  average  over 
all  bands  and  over  all  k  vectors  (i.e.,  for  both  normal  and 
off-axis  incidence).  While  the  specific  heat  and  phonon  phase 
velocity  can  be  determined  quite  easily  from  the  phonon  dis¬ 
persion  curves  of  a  perfect  crystal,  determination  of  the  pho¬ 
non  transmission  coefficient  requires  a  detailed  analysis  of 
the  interface. 

Young  and  Maris40  determined  (t(co))  for  a  nearest- 
neighbor  spring  model  of  a  fee  lattice  with  different  masses 
and  spring  constants  on  the  two  sides;  this  has  recently  been 
extended  to  a  disordered  interface.41  In  both  studies,  it  was 
found  that  ( t(co ))  decreases  almost  uniformly  with  increas- 
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in g  o).  These  results  are  rather  different  from  a  naive  imple¬ 
mentation  of  the  acoustic  mismatch  model,  which  would  pre¬ 
dict  a  weak  dependence  on  frequency  up  to  a  cutoff 
frequency  determined  by  the  maximum  frequency  available 
to  the  softer  material.  A  simple  analytically  tractable  model 
that  captures  the  basic  features  of  the  frequency  dependence 
in  the  phonon  transmission  coefficient  would  be  of  great  use. 

The  experiments  described  in  Sec.  II B  point  to  the  po¬ 
tential  capabilities  of  methods  in  which  phonon  packets  of 
well-defined  wave  vector  and  polarization  are  launched  to¬ 
wards  individual  interfaces.  Recently,  Schelling  et  al.42  per¬ 
formed  atomic-level  simulations  of  phonon-interface  inter¬ 
actions.  The  basic  idea  is  to  launch  phonon  wave  packets  of 
well-defined  polarization  and  wave  vector  towards  interfaces 
and  to  analyze  the  transmitted  and  reflected  phonons.  For 
example,  Fig.  4  shows  four  snapshots  in  the  evolution  of  a 
longitudinal  acoustic  phonon  as  it  approaches,  interacts  with, 
and  then  is  partially  transmitted  through  and  partially  re¬ 
flected  by  an  interface.42  Analysis  of  simulations  such  as  this 
will  allow  the  atomic-level  details  of  phonon-interface  inter¬ 
actions  to  be  elucidated. 

As  discussed  above,  computational  methods  based  on 
applications  of  Fourier’s  Law  and  the  Boltzmann  transport 
equation  are  well  developed.  While  MD  simulation  holds 
promise  of  providing  detailed  atomic  and  spectral  informa¬ 
tion  on  heat  flow  and  phonon  dynamics,  MD  will  not  replace 
BTE  and  Fourier ’s-law  approaches  to  heat  transfer  problems 
in  complex  materials  structures  and  devices.  An  important 
challenge  is  therefore  the  development  of  techniques  to  in¬ 
corporate  the  insights  of  atomic-level  approaches  into  the 
high-level  approaches;  success  in  this  can  be  expected  to 
result  in  BTE  and  Fourier’s  Law  approaches  of  significantly 
higher  fidelity. 


D.  What  is  temperature? 

In  the  above  discussions  of  atomic-level  simulation  we 
have  swept  an  important  conceptual  issue  under  the  rug.  The 
usual  definition  of  temperature  is  related  to  the  average  en¬ 
ergy  of  a  system  of  particles.  This  definition  is  for  a  system 
in  equilibrium,  and  works  even  for  nanoscale  systems.  How¬ 
ever,  our  interest  is  in  the  transport  of  heat  through  nanoscale 
systems.  Can  temperature,  which  is  an  equilibrium  concept, 
still  be  invoked  in  a  nonequilibrium  process  such  as  heat 
flow?  The  answer  is  affirmative  for  macroscopic  systems: 
they  are  so  large  that  one  can  define  a  local  temperature  in 
each  region  in  space.  This  local  temperature  will  vary  from 
region  to  region.  Then  one  finds,  for  example,  that  the  heat 
current  is  proportional  to  —  V7\  The  question  “ What  is  Tem¬ 
perature ?”  is  really  a  question  about  the  size  of  the  regions 
over  which  a  local  temperature  can  be  defined.  In  many 
semiconductor  superlattices,  the  layer  thickness  is  2-5  nm. 
Are  these  regions  large  enough  to  define  a  temperature? 

As  discussed  above,  there  have  been  several  molecular- 
dynamics  simulations  of  heat  flow  through  grain  boundaries. 
The  most  important  conceptual  problem  raised  by  these 
simulations  is  how  to  define  the  temperature  at  different 
planes  in  the  simulation  cell.  The  simulations  calculate  the 
position  R/(X)  and  velocity  \i(tn)  of  each  atom  at  each  time 


step.  It  is  a  relatively  simple  task  to  store  the  velocities,  and 
compute  an  average  kinetic  energy  over  N  steps  in  time 


N 
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It  is  found  that  the  averaging  must  be  done  over  very  long 
times  to  obtain  good  statistical  average  kinetic  energies. 

How  is  this  simple  quantity  converted  to  temperature? 
The  numerical  simulations  show  a  sharp  drop  in  the  average 
kinetic  energy  at  the  grain  boundary,  as  shown  in  Fig.  3. 
Here  each  point  corresponds  to  the  average  kinetic  energy  in 
a  slice  of  a  few  unit-cell  thicknesses.  The  sharp  drop  in  ki¬ 
netic  energy,  at  the  grain  boundary,  is  consistent  with  the 
idea  of  Kapitza  resistance  RK  which  assumes  the  heat  current 
through  the  boundary  is  proportional  to  the  temperature  drop 
Q  =  \T/RK .  In  order  to  find  RK ,  the  kinetic  energy  must  be 
converted  to  a  temperature  scale. 

Because  the  calculations  are  entirely  classical,  with  the 
motion  of  each  silicon  atom  found  from  Newton’s  law,  the 
most  obvious  temperature  scale  is  provided  by  the  classical 
expression 


(4) 


The  temperature  at  each  atom  at  site  R*  is  given  by  T{ .  This 
is  the  procedure  adopted  in  the  MD  simulations  described 
above. 

However,  quantum  mechanics  provides  another  possible 
definition.  The  collective  excitations  of  the  atomic  motions 
are  phonons.  They  have  a  frequency  cox(q)  which  depends 
upon  the  wave  vector  q  and  polarization  X  of  the  phonon. 
The  simulations  of  Fig.  3  were  performed  on  silicon  that  has 
six  polarizations  for  each  wave  vector.  In  quantum  mechan¬ 
ics  the  average  kinetic  energy  of  the  phonons  is 

1 2 mv ^)  =  4N^q  W(q)£W2(Vq) 

X  1  eA^q)/kBTi_l  +  1  j  .  (5) 


where  N  is  the  number  of  unit  cells  which  is  also  the  number 
of  q  points.  The  polarization  vectors  are  .  They  are  nor¬ 
malized  so  that  Ex£(i)2  =  3.  In  the  limit  of  high  temperature, 
where  kBT>ha)x( q)  the  Bose-Einstein  occupation  factors 
are  approximated  as 

2  i  ^  kBT{ 

efta,x(W,_1+l~2^-^y.  (6) 

In  this  high  temperature  limit,  the  quantum  relation  Eq.  (5) 
becomes  identical  to  the  classical  Eq.  (4).  The  problem  with 
many  simulations  is  that  they  are  not  in  the  high  temperature 
limit.  Silicon  has  optical  phonons  of  very  high  energy 
(62  me 750  K).  For  any  temperature  between  300  and 
1000  K,  the  two  definitions  of  temperature  [Eqs.  (4)  and  (5)] 
gave  very  different  values  for  T{ . 

Another  possible  definition  of  temperature  is  to  remove 
the  zero  point  motion  from  the  quantum  definition,  since  it  is 
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not  a  classical  effect.  This  option  now  makes  three  possible 
definitions  of  local  temperature.  Which  is  the  correct  defini¬ 
tion? 

An  important  issue  is  the  size  of  the  region  over  which 
temperature  is  defined.  The  classical  definition  is  entirely 
local,  and  one  can  define  a  temperature  for  each  atom  or 
plane  of  atoms.  For  the  quantum  definition,  the  length  scale 
is  defined  by  the  mean-free-path  €Xq  of  the  phonon.  If  two 
regions  of  space  have  a  different  temperature,  then  they  have 
a  different  distribution  of  phonons.  The  phonons  can  change 
their  distribution  by  scattering.  The  most  important  scattering 
is  the  anharmonic  process  in  which  one  phonon  divides  into 
two,  or  two  combine  to  one.  This  process  occurs  on  the 
length  scale  of  the  mean  free  path.  A  local  region  with  a 
designated  temperature  must  be  larger  than  the  phonon  scat¬ 
tering  distance.  However  low  frequency  phonons  have  a  long 
€Xq,  and  high  frequency  phonons  have  a  short  €Xq.  For  the 
phonons  which  carry  most  of  the  heat,  one  can  plausibly 
define  an  average  mean-free-path:  This  €Xq  is  typically  larger 
than  the  length  of  the  MD  simulation  cell:  this  case  is  called 
the  Casimir43  limit.  One  can  scale  the  results  to  an  infinite 
sized  system  using  the  method  of  Oligschleger  and  Schon.44 

This  phonon  viewpoint  of  temperature  implies  that  tem¬ 
perature  cannot  be  defined  for  a  particular  atom,  or  a  plane  of 
atoms.  In  particular,  there  should  not  be  an  abrupt  variation 
in  temperature  between  a  plane  of  atoms.  Although  this  defi¬ 
nition  seems  quite  reasonable,  it  makes  the  numerical  results 
shown  in  Fig.  3  quite  puzzling.  The  MD  simulations  by  dif¬ 
ferent  groups  do  show  an  abrupt  change  in  the  kinetic  energy 
of  a  plane  of  atoms  at  the  twin  boundary.  Regardless  of 
which  temperature  scale  is  adopted  [Eqs.  (4)  or  (5)],  a  graph 
of  temperature  versus  distance  will  show  an  abrupt  change.  A 
possible  resolution  for  this  puzzle  is  that  a  grain  or  twin 
boundary  may  form  a  natural  boundary  for  a  region  of  tem¬ 
perature.  The  statement  that  temperature  cannot  be  defined 
within  a  scale  of  distance  given  by  €Xq  may  not  apply  across 
grain  boundaries.  The  boundaries  may  provide  natural  limits 
to  the  regions  of  temperature.  Even  if  one  adopts  this  hypoth¬ 
esis,  it  still  means  that  temperature  cannot  vary  within  a 
grain,  or  within  a  superlattice  layer,  on  a  scale  smaller  than 
€Xq.  If  the  layer  thickness  of  the  superlattice  is  less  than 
€Xq,  then  one  cannot  define  T(z)  within  this  layer.  The 
whole  layer  is  probably  at  the  same  temperature.  This  point 
is  emphasized,  since  all  theories  of  heat  transport  in  super¬ 
lattices  have  assumed  that  one  could  define  a  local  tempera¬ 
ture  T(z)  within  each  layer. 

The  results  in  Fig.  3  may  be  measurable.  X-ray  scatter¬ 
ing  at  synchrotron  sources  can  measure  the  lattice  constant  of 
atoms  in  each  atomic  layer.45  The  lattice  constant  can  be 
converted  to  a  local  temperature  scale  for  that  plane.  They 
should  be  able  to  determine  whether  this  quantity  has  an 
abrupt  change  at  an  interface  while  heat  is  flowing.  Although 
this  discussion  of  “What  is  Temperature ?”  has  been  cast  in 
the  framework  of  MD  simulations,  the  issues  are  more  gen¬ 
eral.  One  of  the  major  issues  of  thermal  transport  in  nanos¬ 
cale  systems  is  whether  temperature  can  be  defined  locally.  If 
it  cannot,  then  how  is  transport  calculated? 

Another  approach  is  to  adopt  a  Landauer  formalism,  and 
assume  there  are  two  thermal  reservoirs  at  known  tempera¬ 


tures,  and  consider  the  ballistic  flow  of  heat  between  them. 
This  does  not  work  for  most  nanoscale  devices  since  the 
values  of  €Xq  are  much  smaller  than  the  distance  between  the 
reservoirs.  So  heat  flow  is  diffusive  rather  than  ballistic.  How 
do  you  calculate  diffusive  heat  flow  without  a  local  tempera¬ 
ture? 

As  an  example  of  the  importance  of  temperature,  we  will 
comment  on  two  related  theoretical  calculations.  Stoner  and 
Maris10  measured  the  Kapitza  resistance  between  diamond 
and  other  crystals.  Between  diamond  and  Pb  the  heat  flow 
was  much  higher  than  could  be  explained  by  a  purely  pho¬ 
non  conduction.  Huberman  and  Overhauser46  suggested  that 
the  conduction  electrons  in  Pb  were  playing  an  essential  role 
in  the  boundary  resistance.  Phonons  from  diamond  carry  heat 
to  the  interface,  and  electrons  in  Pb  carry  it  away.  There  have 
been  two  different  calculations  of  this  process,  one  by  Hu¬ 
berman  and  Overhauser,46  and  one  by  Sergeev.47  The  two 
theories  are  very  different  since  the  two  groups  made  differ¬ 
ent  assumptions  regarding  temperature.  Sergeev  assumed 
that  there  was  an  abrupt  change  in  temperature  at  the  bound¬ 
ary  between  diamond  and  Pb.  His  assumption  agrees  with 
the  conventional  view  of  Kapitza  resistance  which  assumes  a 
temperature  step.  Huberman  and  Overhauser  assumed  tem¬ 
perature  could  not  be  defined  locally,  and  there  was  no  step 
in  temperature.  The  interface  region  was  at  one  temperature. 
They  used  a  Landauer  formalism  to  compute  the  energy  ex¬ 
change  between  electrons  and  phonons.  Which  of  these  two 
different  theories  is  better  depends  on  whether  a  grain 
boundary  can  support  an  abrupt  change  in  temperature. 

E.  Transport  theory 

Most  theories  of  transport  in  solids  employ  the  Boltz¬ 
mann  transport  equation  (BTE).  For  both  electron  and  pho¬ 
non  transport,  the  form  of  the  equation,  and  the  form  of  the 
various  scattering  mechanisms,  are  very  well  known.  This 
theory  can  explain,  in  bulk  homogeneous  materials,  the  de¬ 
pendence  of  the  electrical  conductivity,  the  thermal  conduc¬ 
tivity,  and  the  Seebeck  coefficient,  on  temperature,  impurity 
content,  isotope  scattering,  and  quantum  confinement. 

However,  the  Boltzmann  equation  treats  the  electrons 
and  phonons  as  classical  particles.  One  is  solving,  say  for 
phonons,  for  the  density  /[r,ft>x(q),f]  of  excitations  with 
polarization  X,  wave  vector  q,  frequency  o>x(q)  at  point  r  at 
time  t.  The  wave  nature  of  the  excitation  is  neglected,  as  is 
any  interference  phenomena  caused  by  the  wave  nature  of 
the  phonons.  Furthermore,  the  scattering  rates  in  the  BTE  are 
calculated  assuming  the  system  is  only  slightly  perturbed 
from  equilibrium.  The  solutions  to  the  BTE  assume  the  ex¬ 
istence  of  a  local  temperature.  Neither  of  these  two  assump¬ 
tions  may  be  valid  in  nanoscale  devices. 

Wave  interference  becomes  important  in  nanoscale  de¬ 
vices.  The  wavelengths  of  the  phonons  are  similar  to  the 
length  scale  of  the  microstructure.  At  room  temperature,  in 
most  solids,  all  of  the  phonons  states  in  the  Brillouin  zone 
are  involved  in  the  transport.  Their  wavelengths  span  the 
range  from  atomic  dimensions  to  the  size  of  the  sample. 
Most  of  the  actual  heat  is  carried  by  the  phonons  with  wave 
lengths  of  a  few  nanometers.  On  this  scale  wave  interference 
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becomes  important  in  transport.  Similar  considerations  apply 
to  electron  transport  in  nanoscale  devices.  Since  the  tradi¬ 
tional  Boltzmann  equation  cannot  deal  with  these  phenom¬ 
ena,  other  theoretical  approaches  are  required. 

One  method  is  to  calculate  the  phonon  states  for  the 
actual  microstructure,  and  use  those  states  in  the  BTE.  For 
transport  in  a  superlattice,  the  frequencies  o)x(q )  are  those  of 
the  phonons  in  the  superlattice.  This  approach  inserts  in 
some  of  the  wave  phonmena.  The  phonon  states  depend  sig¬ 
nificantly  on  the  interference  due  to  scattering  from  the  mul¬ 
tiple  layer  boundaries  of  the  superlattice.  One  problem  with 
this  approach  is  that  it  ignores  an  important  phenomenon. 
Most  scattering  effects  contribute  to  a  finite  mean-free-path 
which  has  a  major  effect  on  the  phonon  frequencies  in  the 
superlattice.  Solving  this  problem  using  the  BTE  is  also 
complicated  by  phonon  band  folding.  There  are  many  super¬ 
lattice  bands,  phonons  scatter  between  them,  and  the  BTE 
becomes  a  matrix  equation  of  large  dimension.  This  calcula¬ 
tion  is  quite  ambitious  and  has  never  been  done.  This  does 
not  include  all  wave  interference  phenomena. 

Another  possibility  is  to  use  the  quantum  Boltzmann 
equation48  to  solve  for  the  distribution  f(R,t:co,q).  Com¬ 
pared  to  the  BTE,  there  is  one  more  vector  variable  in  the 
argument  that  makes  the  solutions  more  complex.  This  more 
fundamental  equation  does  include  wave  information.  How¬ 
ever,  it  is  difficult  to  solve  and  is  seldom  used.  One  problem 
is  the  choice  of  boundary  conditions  to  apply  on  the  surface 
of  the  system  in  order  to  connect  it  to  a  heat  reservoir.  This 
issue  has  been  much  discussed  in  the  literature  on  electron 
flow  in  nanoscale  devices.49,50 

If  one  throws  out  the  concept  of  a  local  temperature,  so 
that  the  BTE  cannot  be  used,  then  there  is  no  known  way  to 
solve  this  important  problem.  Many  semiconductor  devices 
are  dependent  upon  the  electrical  currents  provided  by  the 
electrons.  The  heat  currents  have  components  from  both 
electrons  and  phonons.  Presently  there  is  no  accurate  way  to 
model  the  heat  flows  in  these  nanoscale  systems  while  in¬ 
cluding  the  exchange  of  heat  between  electrons  and  phonons. 
A  major  theoretical  challenge  is  to  invent  a  new  method  of 
solving  this  and  related  transport  problems. 

III.  PHONONS  AND  THERMAL  TRANSPORT  IN 
NANOSCALE  STRUCTURES  AND  DEVICES 

A.  Silicon  films  and  devices 

1.  Conduction  along  semiconducting  monolayers 

Phonon  conduction  along  single-  and  poly  crystalline 
monolayers  is  important  for  silicon-on-insulator  (SOI) 
circuits,51  novel  3D  multilayer  electronics,52  and  MEMS 
with  semiconducting  or  dielectric  membranes.53  Silicon 
monolayers  have  received  the  most  attention,  either  as  sus¬ 
pended  membrane  (a  geometry  relevant  for  many  thermal 
microsensors)  or  embedded  between  amorphous  silicon  di¬ 
oxide  films  (as  in  SOI  transistors).  In  the  second  geometry, 
the  surrounding  silicon  dioxide  films  act  as  diffuse  phonon 
emitter  absorbers  and  provide  a  boundary  condition  that  re¬ 
sembles  phonon  scattering  on  the  boundaries  of  a  membrane 
with  roughness  above  a  few  angstroms.  This  situation  differs 
from  that  in  superlattices  with  comparable  phonon  mean  free 
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FIG.  5.  Overview  of  thermal  conductivity  data  for  silicon  films  with  varying 
microstructural  quality  and  impurity  concentrations  after  Ref.  58.  The  data 
illustrate  the  relative  importance  of  phonon  scattering  on  film  and  grain 
boundaries  and  impurities.  The  data  for  undoped  single-crystal  bulk  samples 
are  taken  from  Ref.  59.  Both  single  crystal  films  have  thickness  of  3  jmm, 
and  for  the  doped  film,  a  boron  concentration  of  1019  cm-3.  (Ref.  55).  The 
undoped  and  doped  (1.6X  1019  cm-3  boron)  polycrystalline  films  have  grain 
sizes  of  200  and  350  nm,  respectively  (Refs.  57  and  58). 


path  on  both  sides  of  each  interface,  for  which  conduction  is 
influenced  by  the  coupled  ballistic  phonon  transport  in  adja¬ 
cent  layers,  see  Sec.  IV.  (A  macroscopic  version  of  this  ge¬ 
ometry  has  been  used  by  Pohl  and  co-workers54  to  study 
phonon  scattering  in  thin  films  at  T<  1  K. ) 

Figure  5  shows  that  for  single-crystal  silicon  films, 
phonon-interface  scattering  strongly  reduces  the  effective 
thermal  conductivity  at  low  temperatures.  For  single-crystal 
silicon  layers  doped  at  concentrations  higher  than  1017  cm-3 
the  conductivity  is  reduced  by  scattering  on  impurities  and 
free  electrons.55  Figure  5  shows  that  this  is  particularly  im¬ 
portant  at  low  temperatures.  Measurements  at  20  K  on  layers 
of  thickness  3  jmm  with  phosphorus  and  boron  concentrations 
of  1.0X  1018  cm-3  yielded  reductions  by  approximately  two 
and  four,  respectively,  with  the  difference  resulting  from  the 
disparity  in  mass  of  the  two  impurity  types.  At  room  tem¬ 
perature,  measurements  along  silicon  films  of  thickness 
down  to  74  nm  showed  a  reduction  of  the  conductivity  by  a 
factor  of  two  compared  to  bulk  data. 

In  poly  crystalline  silicon  films,  phonon  scattering  on 
grain  boundaries  and  related  defects  dominates  over  scatter¬ 
ing  from  the  surface  or  the  film/substrate  interface.  Poly  sili¬ 
con  conductivities  are  reduced  compared  to  those  of  pure 
crystalline  films,56  but  the  roles  of  impurities  and  grains  are 
difficult  to  separate.  Recent  work  studied  films  with  grain 
sizes  between  300  and  500  nm,  extracted  using  transmission 
electron  microscopy,  and  boron  or  phosphorus  concentra¬ 
tions  up  to  4.1X1019  cm-3,  measured  using  secondary  ion 
mass  spectroscopy  (SIMS).57,58  Figure  5  shows  that  poly  sili¬ 
con  conductivities  are  reduced  strongly  at  all  temperatures 
compared  to  similarly  doped  single-crystal  silicon  layers, 
which  illustrates  the  importance  of  grain  boundary  scatter¬ 
ing.  The  grain  structure  and  spatial  impurity  distribution  de¬ 
pend  strongly  on  the  temperature  and  duration  of  annealing. 
Annealing  of  the  doped  polysilicon  film  and  the  resulting 
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FIG.  6.  Experimental  data  for  the  effective  room-temperature  thermal  con¬ 
ductivity  along  silicon  films  of  thickness  near  100  nm  (Ref.  64).  The  data 
are  compared  with  predictions  based  on  the  Boltzmann  transport  equation 
engaging  different  assumptions  about  the  frequency  dependence  of  phonon 
relaxation  times  and  the  relative  contributions  of  the  phonon  branches  to 
conduction. 


increase  in  grain  size  yields  a  higher  thermal  conductivity 
than  that  of  the  undoped  film. 

The  conduction  size  effect  along  crystalline  and  poly¬ 
crystalline  monolayers  at  low  temperatures  can  be  effectively 
predicted  using  a  semiclassical  approach,  based  on  solutions 
to  the  Boltzmann  transport  equation  for  conduction  along 
films60  and  based  on  a  spectrally-resolved  form  of  the  ex¬ 
pression  k=  Cu€/3.61  In  this  expression,  C  is  the  heat  ca¬ 
pacity  per  unit  volume  and  v  and  €  are  the  phonon  velocity 
and  mean  free  path,  respectively.  Polysilicon  data  are  effec¬ 
tively  predicted  below  about  10  K  by  reducing  the  mean  free 
path  using  Matthiessen’s  rule62  and  a  grain-boundary  free 
path  given  by  dGIB,  where  dG  is  the  grain  size  and  B  is  a 
dimensionless  constant  governed  by  the  grain  shape  and 
grain  boundary  reflection  coefficient.  Polycrystalline  dia¬ 
mond  films  are  worthy  of  brief  discussion  here  because  they 
illustrate  another  way  in  which  grain  boundaries  can  reduce 
phonon  conduction.  The  free  path  dG/B  is  overly  simplistic 
for  diamond  film  data  owing  to  the  large  concentration  of 
other  imperfections.  Because  imperfections  populate  prima¬ 
rily  at  grain  boundaries,  the  phonon  scattering  rate  can  still 
be  coupled  to  the  grain  size  dG  using  a  semiclassical 
approach.63 

The  anisotropy  and  nonlinearity  of  phonon  dispersion 
relationships  complicate  predictions  above  100  K  for  GaAs, 
Ge,  and  AlAs  films,  and  above  150  K  for  silicon  films  owing 
to  its  higher  Debye  temperature.  The  predictions  become 
more  challenging  because  phonon  dispersion  induces  differ¬ 
ences  in  transport  among  the  various  phonon  modes  and  di¬ 
rections.  Figure  6  illustrates  this  issue  for  transport  along 
silicon  films  of  thickness  near  and  below  100  nm.  Predic¬ 
tions  are  performed  under  a  greybody  assumption,  which 
overestimates  the  fraction  of  phonons  contributing  to  heat 
transfer  and  therefore  underestimates  the  phonon  mean  free 
path  and  the  thermal  conductivity  size  effect  for  a  film  of 
given  thickness.  Predictions  are  also  provided  with  spectrally 
varying  mean  free  paths  under  the  competing  assumptions 
about  which  phonon  modes  are  dominating  conduction  in  the 
material. 


2.  Monte  Carlo  solutions  of  the  Boltzmann  transport 
equation 

The  problem  of  phonon  transport  considering  dispersion 
effects  has  been  addressed  using  Monte  Carlo  simulations, 
which  are  essentially  an  alternative  approach  to  solving  the 
Boltzmann  transport  equation.  The  Boltzmann  transport 
equation  (BTE)  describes  heat  flow  by  phonons  in  nonmetal- 
lic  solids,  when  wave  effects  are  negligible.  The  BTE,  in  its 
most  general  form,  is  difficult  to  solve  analytically  for  real¬ 
istic  phonon  dispersion  and  density  of  states  as  well  as  for 
transitions  between  phonon  polarizations  and  for  irregular 
geometries.65  Recent  work-by  Majumdar,26  Chen  and  Tien,66 
Goodson,63  and  Chen67  have  used  analogies  from  radiative 
heat  transfer  and  presented  solution  strategies  for  both  the 
diffusion  and  the  ballistic  limit.  In  such  calculations,  how¬ 
ever,  phonon  dispersion  was  neglected,  and  a  single  “aver¬ 
age”  polarization  branch  was  considered.  Furthermore,  such 
calculation  strategies  can  only  be  adopted  for  considerably 
simplified  geometries.  If  all  these  assumptions  are  removed, 
it  is  extremely  difficult  to  solve  the  BTE  for  phonons  by  a 
deterministic  approach  in  an  arbitrary  geometry.  This  is  be¬ 
cause  the  number  of  independent  variables  is  too  large,  and 
would  render  any  kind  of  discretization  scheme  too  complex 
to  be  practical.  Furthermore,  the  phonon-phonon  scattering 
events  are  difficult  to  incorporate  without  a  simplifying  re¬ 
laxation  time  approximation. 

One  approach  to  capture  the  complexity  of  phonon  trans¬ 
port  in  the  particle  regime  is  to  use  Monte  Carlo  or  stochastic 
simulations  to  solve  the  BTE.  Over  the  last  decade  there  has 
been  tremendous  advancement  in  the  development  of  Monte 
Carlo  solution  techniques  for  the  BTE  for  electrons  and  holes 
in  semiconductors.68-72  However,  there  have  been  very  few 
reports  of  using  Monte  Carlo  simulations  for  phonon 
transport.73,74  These  studies  have  either  assumed  the  Debye 
approximation  for  density  of  states  or  did  not  consider 
phonon-phonon  interactions.  More  recently,  Mazumder  and 
Majumdar75  considered  phonon  dispersion  as  well  as  various 
phonon  scattering  mechanisms  to  study  heat  transport  in 
complex  geometries  and  to  predict  the  thermal  conductivi¬ 
ties.  The  result  showed  that  by  fitting  one  parameter  using 
experimental  data  of  thermal  conductivity  at  one  tempera¬ 
ture,  predictions  of  the  thermal  conductivity  of  silicon  agreed 
well  with  experimental  data  over  a  wide  temperature  range. 
In  addition,  it  was  able  to  capture  phonon  transport  in  both 
the  ballistic  and  diffusive  limits.  Despite  the  success  of  this 
approach,  one  must  be  aware  of  the  challenges  in  applying 
this  approach  more  generally  to  problems  in  nanoscale  ther¬ 
mal  transport. 

First,  it  must  be  recognized  that  it  is  difficult  to  capture 
all  the  selection  rules  for  transitions  between  polarizations 
during  phonon-phonon  interactions.  For  example,  the  tran¬ 
sitions  between  longitudinal  and  transverse  acoustic  phonons 
as  well  as  those  between  acoustic  and  optical  phonons  are 
very  important  in  energy  dissipation  and  thermal  transport. 
Mazumder  and  Majumdar75  ignored  the  selection  rules  and 
obtained  local  phonon  distributions  in  different  polarizations 
based  on  a  quasiequilibrium  condition.  For  highly  nonequi¬ 
librium  cases,  this  approach  cannot  be  applied.  Hence,  it  is 
extremely  important  to  develop  schemes  for  capturing  tran- 


Downloaded  18  Nov  2008  to  171.64.49.29.  Redistribution  subject  to  AIP  license  or  copyright;  see  http://jap.aip.org/jap/copyright.jsp 


J.  Appl.  Phys.,  Vol.  93,  No.  2,  15  January  2003 


Appl.  Phys.  Rev.:  Cahill  et  al. 


803 


sitions  between  polarizations  during  phonon-phonon  inter¬ 
actions.  Second,  the  relaxation  times  for  phonon-phonon  in¬ 
teractions  have,  in  the  past,  been  used  only  for  bulk  solids 
and  only  for  the  most  common  materials.  These  were  ob¬ 
tained  from  approximate  analytical  solutions  that  require  fit¬ 
ting  parameters.  For  example,  the  phonon  scattering  relax¬ 
ation  times  for  silicon  are  generally  obtained  from 
Holland.61,76  These  cannot  be  used  when  silicon  is  nanostruc- 
tured  which  can  modify  the  phonon  dispersion  relation.  For 
many  other  materials,  relaxation  times  are  not  available  even 
for  bulk  solids.  Hence,  there  is  a  need  to  develop  approaches 
of  predicting  the  relaxation  times  for  arbitrary  phonon  dis¬ 
persion  relations.  This  would  allow  one  to  study  energy  dis¬ 
sipation  and  thermal  transport  in  nano  structured  solids  such 
as  quantum  dots,  nano  wires,  films,  and  superlattices,  where 
the  phonon  dispersion  relations  are  modified  due  to  confine¬ 
ment  or  wave  interference.  If  both  these  challenges  are  over¬ 
come,  Monte  Carlo  simulations  of  phonon  transitions  and 
transport  can  be  an  extremely  effective  tool  to  study  nanos¬ 
cale  energy  dissipation  and  thermal  transport. 

3.  Nanoscale  hotspots  in  devices 

The  phonon-boundary  scattering  discussed  above  for 
silicon  monolayers  is  only  the  first  among  three  nanoscale 
effects  which  complicate  heat  conduction  from  metal-oxide- 
semiconductor  (MOS)  transistors.  Two  additional  nanoscale 
effects  result  from  the  small  size  of  the  region  of  electron 
energy  transfer  to  the  lattice.  Before  describing  these  effects, 
it  is  important  to  note  that  all  three  are  augmenting  the  de¬ 
parture  from  equilibrium  of  the  phonon  distribution  func¬ 
tions,  rendering  the  usefulness  of  a  phonon  “temperature” 
more  questionable.  Theoretical  treatments  of  these  effects 
use  temperatures  either  for  the  phonon  system  as  a  whole  or 
for  fractions  of  the  phonons  sharing  comparable  group  ve¬ 
locities;  this  approach  can  at  best  be  viewed  as  a  shorthand 
placeholder  for  energy  density.  This  situation  may  eventually 
be  remedied  through  detailed  phonon  transport  analysis  or 
molecular  dynamics,  although  the  computational  demands 
will  be  excessive  owing  to  the  three-dimensional  nature  of 
the  problem. 

Electrons  traveling  between  the  source  and  drain  are  ac¬ 
celerated  by  the  electric  fields  and  generate  primarily  optical 
phonons.  Simulations  of  electron  transport  using  drift  diffu¬ 
sion  and  the  more  rigorous  Monte  Carlo  or  hydrodynamic 
approaches77,78  illustrate  that  rates  of  phonon  generation  are 
sharply  peaked  on  the  drain  side.  The  region  with  strong 
lattice  heat  generation  has  dimensions  near  tens  of  nanom¬ 
eters  and  is  much  smaller  than  the  channel  length  (50-150 
nm  depending  on  the  technology  generation  of  interest)  and 
the  phonon  mean  free  path,  €^250  nm  according  to  Ref.  64. 
This  second  nanoscale  effect,  which  we  will  call  hotspot  bal¬ 
listic  emission  effect,  causes  the  temperature  of  the  phonon 
“hotspot”  to  be  larger  than  that  predicted  using  diffusion 
theory,  by  a  factor  comparable  with  €/r  where  r  is  the 
hotspot  dimension.79  While  the  precise  distribution  of  the 
generated  phonons  between  the  optical  and  acoustic  modes 
depends  strongly  on  the  energies  of  scattering  electrons,  it  is 
clear  that  the  generated  phonons  have  group  velocities  far 
below  the  velocity  characteristic  of  heat  conduction  by  the 


FIG.  7.  Experimental  data  for  the  phonon  temperature  rise  within  a  resistive 
heater  in  silicon  with  dimensions  (at  low  temperatures)  smaller  than  the 
phonon  mean  free  path.  The  data  are  compared  with  predictions  based  on 
diffusion  theory,  the  two-fluid  phonon  Boltzmann  transport  equation,  and  a 
simple  resistance  model  based  on  the  ratio  of  the  square  of  the  mean  free 
path  and  the  effective  cross-sectional  area  of  the  heat  source  (Ref.  81). 


lattice.  The  generation  of  high  concentrations  of  optical 
phonons  within  a  region  with  dimensions  small  compared  to 
the  mean  free  path  causes  a  third  nanoscale  effect,  which  in 
the  present  manuscript  we  will  call  the  hotspot  far-from- 
equilibrium  effect. 

Sverdrup80,81  examined  this  third  nanoscale  effect  ex¬ 
perimentally  using  a  dedicated  micro  structure  and  theoreti¬ 
cally  by  breaking  the  phonons  into  two  propagating  and  sta¬ 
tionary  groups.  The  theoretical  treatment  solved  the 
Boltzmann  transport  equation  for  the  propagating  phonons 
together  with  a  reservoir  energy- storage  equation  for  the  sta¬ 
tionary  phonons.  Stationary  phonons  were  generated  by  elec¬ 
tron  scattering,  with  a  spatial  distribution  taken  from  an  in¬ 
dependent  analysis  of  the  electron  and  hole  transport 
problem  in  the  transistor.  Simulations  of  a  SOI  transistor 
with  channel  length  150  nm  concluded  that  the  peak  lattice 
temperature  (a  weighted  average  of  the  temperatures  of  the 
stationary  and  propagating  phonons)  was  elevated  by  more 
than  a  factor  of  two  compared  to  predictions  based  only  on 
the  heat  diffusion  equation.  This  difference  resulted  in  part 
from  the  localized  transfer  of  energy  from  phonons  with  rela¬ 
tively  small  group  velocities  to  those  at  lower  energies, 
which  dominate  heat  conduction. 

This  third  nanoscale  effect  was  examined  experimentally 
using  localized  heating  by  a  doped  silicon  bridge  within  a 
suspended  silicon  membrane.  Experimental  data  and  predic¬ 
tions  for  the  temperature  of  the  bridge  are  provided  in  Fig.  7. 
The  cross-sectional  dimensions  of  the  doped  heater  were 
~  1  /mm  and  much  larger  than  the  phonon  mean  free  path  at 
room  temperature,  such  that  all  of  the  predictions  and  the 
data  are  consistent  with  simple  diffusion  theory.  At  lower 
temperatures,  where  the  phonon  mean  free  path  is  longer,  the 
dramatic  departure  from  equilibrium  induced  by  optical  pho¬ 
non  generation  within  the  hotspot  (far-from-equilibrium  ef¬ 
fect)  becomes  important.  Predictions  at  low  temperatures  us¬ 
ing  the  Boltzmann  equation  with  the  stationary  and 
propagating  phonon  groups  agree  reasonably  well  with  the 
data,  and  are  consistent  with  a  simple  closed-form  expression 
in  which  the  excess  temperature  rise  scales  with  (€/r)2. 
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The  hotspot  effect,  discussed  above  for  silicon,  could 
very  well  be  much  larger  in  optoelectronic  devices  because 
of  the  phonon  band  structure  in  the  associated  semiconduct¬ 
ing  materials.  In  polar  semiconductors,  such  as  GaAs,  the 
carriers  interact  predominantly  with  LO  phonons  of  wave 
vector  0.82  Thus,  carrier  cooling  leads  to  the  generation 
of  nonequilibrium  modes  for  which  the  group  velocity  is 
exceedingly  small,  i.e.,  to  a  nanoscale  hot  spot.  This  so- 
called  hot-phonon  effect  is  believed  to  be  responsible  for  the 
reduction  of  the  cooling  rate  of  electrons  at  densities  higher 
than  1017  cm-3.83  It  is  well  known  that  the  presence  of  in¬ 
terfaces  changes  the  character  of  the  optical  phonons  in  that 
they  are  no  longer  purely  LO  or  transverse-optical  (TO)  but  a 
mixture  involving  LO,  TO,  and  interface  modes.84  Since 
mixing  results  in  larger  dispersion,  this  suggests  that  struc¬ 
ture  engineering  can  be  applied  to  attain  large  group  veloci¬ 
ties  to  facilitate  the  removal  of  LO  heat  from  the  device 
active  region.  In  the  bulk,  heat  diffusion  relies  on  the  spon¬ 
taneous  decay  of  the  optical  phonons  into  the  more  movable 
acoustic  phonons.85  Recent  progress  in  the  generation  of  co¬ 
herent  high-frequency  acoustic  modes,86  which  can  be  used 
to  stimulate  LO  decay,  hold  promise  for  applications  in 
nanoscale  heat  removal. 


B.  Nanostructures  with  tailored  phonon  transport 
properties 

Unlike  bulk  materials,  the  phonon  properties  of  semicon¬ 
ductor  nanostructures  and,  in  particular,  the  phonon  fre¬ 
quency,  group  velocity,  spectral  density  as  well  as  the 
strength  of  the  interaction  with  carriers,  can  be  widely  modi¬ 
fied  to  improve  the  performance  of  optoelectronic  and  other 
devices.  In  some  cases,  it  is  desirable  to  increase  the 
electron-phonon  interaction;  an  example  is  the  interband 
semiconductor  laser,  where  fast  relaxation  rates  reduce  hot 
electron  effects.87  On  the  opposite  side,  we  have  devices 
such  as  mid-infrared  emitters  and  detectors  based  on  inter¬ 
subband  transitions,  for  which  a  reduction  of  the  electron- 
phonon  coupling  constant  is  advantageous  because  it  leads  to 
larger  lifetimes  of  the  relevant  excited  states.87  An  important 
example  of  the  effect  of  quantum  confinement  on  scattering 
rates  is  that  of  the  phonon  bottleneck  in  quantum  dots.  This 
refers  to  the  strong  suppression  of  the  electron  relaxation  rate 
which  arises  when  the  energy  spacing  of  the  confined  carri¬ 
ers  exceeds  the  energy  of  the  longitudinal  optical  (LO) 
phonon.88 


C.  Phonon  transport  in  ID  nanostructures 

One-dimensional  (ID)  nanostructures  such  as  carbon 
nanotubes  and  semiconductor  nanowires  have  recently  re¬ 
ceived  a  lot  of  attention.  While  most  of  the  current  research 
is  being  focused  on  electronic  and  optical  properties,  thermal 
transport  is  also  of  interest  for  basic  science  as  well  as  for 
technological  applications.  For  example,  Schwab  et  al ,89 
demonstrated  quantum  thermal  conductance  in  nanofabri- 
cated  ID  nanostructures,  where  the  quantum  of  phonon  con¬ 
ductance,  Gph ,  is  given 


FIG.  8.  A  large  scale  scanning  electron  microscopy  image  of  a  microfabri- 
cated  device  (left).  Two  independent  islands  are  suspended  by  three  sets  of 
250  /mm  long  silicon  nitride  legs  with  Pt  lines  that  connect  the  microther¬ 
mometer  on  the  islands  to  the  bonding  pads,  (right)  Enlarged  image  of 
another  design  of  suspended  islands  with  the  Pt  thermometer  and  heater 
fabricated  by  electron  beam  lithography.  Also  shown  is  an  individual 
multiwall  carbon  nanotube,  14  nm  in  diameter,  bridging  the  two  islands 
(Ref.  94). 


TT2k2BT 

G-*-— •  <7> 

where  kB  is  the  Boltzmann  constant,  T  is  the  absolute  tem¬ 
perature,  and  h  is  the  Planck  constant. 

In  this  regime,  a  ID  nanostructure  behaves  essentially 
like  a  phonon  waveguide  similar  to  optical  ones  for  light. 
Berber90  and  Osman  and  Srivastava91  have  theoretically  pre¬ 
dicted  very  high  thermal  conductivity  for  carbon  nanotubes. 
In  addition,  Hicks  and  Dresselhaus92  have  suggested  that 
quantum  confinement  of  electrons  in  nanowires  could  be 
used  to  manipulate  thermoelectric  properties,  which  could  be 
useful  for  solid-state  energy  conversion  devices.  Hence,  it  is 
important  to  study  thermal  transport  and  energy  conversion 
in  ID  nanostructures. 

On  the  experimental  side,  there  have  been  some  recent 
attempts93  to  measure  the  thermal  conductivity  of  a  collec¬ 
tion  or  a  mat  of  carbon  nanotubes.  However,  due  to  weak 
coupling  between  carbon  nanotubes,  the  measured  thermal 
conductivity  was  found  to  be  much  lower  than  that 
predicted.90  Kim94  developed  a  microdevice  (see  Fig.  8)  con¬ 
taining  two  adjacent  silicon  nitride  (SiNx)  membranes  or  is¬ 
lands  suspended  with  long  thermal  isolation  legs.  Each  island 
contains  a  Pt  resistor  which  can  be  used  as  both  a  heater  and 
a  thermometer.  A  nanowire  can  be  placed  in  between  and 
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FIG.  9.  Measured  thermal  conductivity  of  a  14  nm  diameter  multiwall  car¬ 
bon  nanotube  (MWCN)  (solid  circle),  an  80  nm  diameter  MWCN  bundle 
(solid  triangle),  and  a  200  nm  diameter  MWCN  bundle  (solid  square)  (Ref. 
94).  Data  for  two  vapor-grown  graphite  fibers  (Ref.  95),  one  heat  treated  to 
3000  °C  (open  triangle)  and  one  without  heat  treatment  (open  circle)  are 
included  for  comparison.  The  lines  represent  the  calculated  (Ref.  95)  basal- 
plane  thermal  conductivity  of  graphite,  assuming  temperature-independent 
low-temperature  phonon  mean  free  path  €  =  2.9  jmm  (upper  line)  and  € 
=  3.9  nm  (lower  line). 


bridge  the  two  suspended  islands.  One  of  the  islands  is 
heated  to  a  temperature  Th  ,  which  is  determined  by  measur¬ 
ing  the  resistance  of  the  Pt  wire  on  the  island.  Conduction 
through  the  nanowire  heats  the  sensing  island  to  a  tempera¬ 
ture  Ts ,  which  is  also  measured  by  resistance  thermometry. 
By  noting  the  power  dissipated  by  the  heater,  Qh ,  and  esti¬ 
mating  the  thermal  conductance  of  the  suspended  legs  con¬ 
nected  the  islands  in  the  absence  of  a  nanowire,  the  conduc¬ 
tance  of  the  nanowire  Gn  can  be  calculated. 

Figure  9  shows  the  plot  of  thermal  conductivity  of  a 
multiwall  carbon  nanotube  (MWCN)  as  a  function  of  tem¬ 
perature  measured  by  this  technique.  The  T 2  temperature  de¬ 
pendence  suggests  that  MWCN  behaves  thermally  as  a  2D 
solid;  the  MWCN  has  a  thermal  conductivity  of  about 
3000  W  m_  1  K- 1  at  room  temperature,  which  is  in  close 
agreement  with  the  predictions  of  Osman  and  Srivastava.91 
Note  that  the  measured  value  includes  the  contact  resistance 
between  the  nanotube  and  the  two  islands,  although  it  is 
predicted  to  be  lower  than  the  resistance  of  the  tube  itself. 
Estimation  of  the  contact  resistance  has  remained  a  chal¬ 
lenge.  Nevertheless,  it  is  envisioned  that  this  technique  can 
be  used  to  measure  thermal  conductivities  of  a  variety  of 
semiconducting  or  metallic  nanowires  that  are  technologi¬ 
cally  important. 

Phonon  transport  in  ID  nanostructures  offers  the  oppor¬ 
tunity  to  understand  the  basic  science  of  phonon  dynamics 
and  transport,  while  allowing  the  ability  to  manipulate  ther¬ 
mal  properties.  Phonons  in  nanowires  can  be  different  from 
those  in  bulk  semiconductors  mainly  because  the  dispersion 
relation  could  be  significantly  modified  due  to  confinement 
in  two  directions.  In  addition,  the  presence  of  a  surface  can 
introduce  surface  phonon  modes.  These  result  in  many  dif¬ 
ferent  phonon  polarizations  other  than  the  two  transverse  and 


150  nm -  40  nm  - 


FIG.  10.  Images  of  divots  written  in  a  polymer  film  using  heat  and  force 
applied  by  a  silicon  atomic-force-microscope  (AFM)  cantilever.  The  images 
are  recorded  using  the  same  cantilever  used  during  the  writing  process  and 
a  thermal  reading  technique  involving  heat  conduction  from  the  cantilever 
into  the  substrate  (Ref.  53). 


one  longitudinal  acoustic  branches  found  in  bulk  semicon¬ 
ductors.  Such  changes  in  the  dispersion  relation  can  modify 
the  group  velocity  and  the  density  of  states  of  each  branch. 
The  phonon  lifetime  also  changes  and  this  arises  from  two 
sources.  First,  the  phonon-phonon  interactions  can  change 
because  selection  rules  based  on  energy  conservation  and 
wave-vector  relations  depend  on  the  dispersion  relation.  Sec¬ 
ond,  boundary  scattering  can  be  much  stronger  in  nanowires 
(5-50  nm  diameter)  than  in  bulk  semiconductors.  While 
there  have  been  some  theoretical  studies  on  phonon  conduc¬ 
tion  in  nano  wires, 96,97  there  are  still  several  open  questions 
regarding  phonon  transport  in  ID  nanostructures:  What  are 
the  phonon-phonon  relaxation  times  and  the  selection  rules? 
How  does  the  phonon  gas  couple  to  the  electronic  one?  What 
is  the  role  of  individual  defects  in  a  nanowire  on  phonon 
transport  and  does  the  presence  or  absence  of  single  defects 
alter  overall  thermal  transport?  How  can  nanowires  be  de¬ 
signed  for  very  high  or  low  thermal  conductivities?  It  is  en¬ 
visioned  that  these  questions  can  only  be  answered  by  a 
combination  of  experimental  studies  of  single  nanowires  and 
theoretical  and  computational  (molecular  dynamics  and 
Monte  Carlo  simulation)  studies  of  phonon  transport  in  ID 
nanostructures. 


D.  Nanoscale  3D  conduction  and  mass  transport  in 
polymers 

Section  III  focused  on  heat  transport  in  crystalline  or 
nearly  crystalline  materials.  Recent  developments  in  data 
storage  technology  highlight  important  shortfalls  in  our  un¬ 
derstanding  of  nanoscale  heat  and  mass  transport  in  disor¬ 
dered  materials  such  as  polymers.  In  thermomechanical  data 
storage,  a  heated  AFM  cantilever  forms  a  divot  with  radius 
of  curvature  below  20  nm  in  a  polymer  film  of  thickness  near 
30  nm  on  a  silicon  substrate.53  Figure  10  illustrates  divots 
written  in  a  polymer  film,  as  observed  using  an  AFM-based 
thermal  readback  mechanism.  The  transport  of  heat  and  mass 
in  the  vicinity  of  the  cantilever  tip  governs  the  bit  writing 
rate  and  bit  size;  therefore,  nanoscale  transport  properties  are 
the  critical  figures  of  merit.  Fundamental  study  is  motivated 
by  the  fact  that  the  thickness  of  the  polymer  film  and  the 
radius  of  curvature  of  the  divot  are  more  than  an  order  of 
magnitude  smaller  than  the  radius  of  gyration  of  the  polymer 
in  bulk  form.  Therefore,  heat  flow  and  deformation  within 
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the  polymer  will  be  strongly  influenced  by  the  partial  orien¬ 
tation  of  molecular  strands  in  the  plane  of  the  film. 

While  the  theory  of  polymer  motion  has  evolved  much 
in  the  past  two  decades,98,99  there  is  little  theory  and  data  of 
direct  relevance  to  the  three-dimensional  motion  observed 
during  bit  formation.  The  three-dimensional,  multiscale  (10- 
100  nm)  nature  of  the  conduction  problem  within  the  canti¬ 
lever  and  into  the  polymer  film  may  well  render  it  beyond 
the  reach  of  atomistic  simulations  for  a  number  of  years. 
Furthermore,  semiclassical  phonon  simulations  such  as  the 
Boltzmann  equation  and  Monte  Carlo  techniques,  while  ap¬ 
propriate  in  the  silicon  cantilever  tip,  are  of  little  use  in  the 
disordered  polymer.  Modifications  of  continuum  theory  to 
account  for  anisotropy  and  partial  orientation  of  molecular 
strands  provide  some  help  with  the  heat  transport  problem,100 
and  may  become  more  accurate  when  coupled  with  atomistic 
simulations  of  individual  molecules. 

IV.  THERMAL  TRANSPORT  IN  NANOSTRUCTURED 
MATERIALS 

In  this  section  we  consider  bulk  or  thin  film  materials 
that  contain  a  collection  of  nano  structured  units.  This  in¬ 
cludes  amorphous  or  glassy  materials,  materials  with  nanos¬ 
cale  grains,  as  well  as  nanoporous  materials  that  contain 
voids  on  the  order  of  a  few  nanometers.  Also  included  in  this 
category  are  multilayer  films  as  well  as  epitaxial  superlat¬ 
tices,  which  include  a  collection  of  nanometer-scale  films 
stacked  on  each  other.  The  common  feature  in  these  materi¬ 
als  is  that  one  can  identify  a  nano  structured  unit  as  the  build¬ 
ing  block.  For  example,  in  amorphous  materials  the  local 
ordering  is  on  the  scale  of  lattice  constants  while  for  epitax¬ 
ial  superlattices  the  unit  is  a  superlattice  period. 

While  phonon  dynamics  and  transport  in  a  single  unit  is 
important  and  must  be  studied,  the  collective  behavior  of  a 
large  number  of  units  can  add  to  more  complexity,  i.e.,  there 
may  be  some  collective  modes  which  are  not  found  when 
considering  only  single  units.  For  example,  when  the  phonon 
coherence  length  scales  are  larger  than  the  size  of  the  unit, 
phonon  interference  effects  lead  to  modified  dispersion  rela¬ 
tions.  Hence,  the  appearance  of  both  wave  effects  at  nanos¬ 
cales  and  diffusive  heat  flow  at  bulk  scales  poses  challenges 
in  predicting  thermal  transport  in  these  materials. 

Such  predictions  are  important  because  nanostructured 
materials  find  many  applications.  For  example,  with  the  size 
of  electronic  devices  shrinking  it  is  necessary  to  use  materi¬ 
als  with  low  dielectric  constants.  Amorphous  or  nanoporous 
materials  are  generally  chosen  for  this  purpose,  although 
their  thermal  conductivities  are  very  low  which  results  in 
poor  thermal  management.  Can  materials  with  low  dielectric 
constants  and  high  thermal  conductivity  be  designed?  This 
can  be  answered  if  we  have  a  fundamental  understanding  of 
phonon  dynamics  and  transport  in  such  materials.  Refractory 
ceramic  films  are  also  used  in  thermal  barrier  coatings  in  gas 
turbine  engines.  The  performance  of  such  engines  can  be 
significantly  improved  if  the  thermal  conductivity  of  these 
films  can  be  reduced,  thus  allowing  operation  at  higher  tem¬ 
peratures.  Finally,  epitaxial  superlattices  of  semiconductor 
films  are  under  intense  investigation  for  thermoelectric  en¬ 
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ergy  conversion  devices  because  they  offer  the  ability  to  the 
control  both  electronic  band  structure  and  phonon  transport. 

A.  Amorphous  and  nanoporous  materials 

Heat  transport  is  an  important  problem  to  the  semicon¬ 
ductor  industry,  which  is  manifested  in  at  least  two  key  ar¬ 
eas:  advanced  packaging  and  device  reliability.  In  the  former 
case,  the  challenge  is  how  to  maintain  acceptable  operating 
temperatures  with  an  increasingly  dense  circuit  layout  oper¬ 
ating  at  increasingly  higher  frequencies.  In  short,  the  pack¬ 
age  cooling  requirement  for  silicon  logic  devices  is  projected 
to  be  greater  than  175  W  after  2010.101  The  package  thermal 
problem  is  amplified  when  one  considers  that  the  increas¬ 
ingly  widespread  proliferation  of  logic  chips  requires  flexible 
package  platforms  for  applications  ranging  from  appliances 
and  personal  digital  assistants  at  one  extreme  to  the  high 
compute  node  density  required  in  modern  server  farms.  Hy¬ 
brid,  chip-on-chip,  and  three  dimensional  chip  construction 
adds  to  the  package  complexity  that  must  be  addressed.  Fur¬ 
thermore,  as  performance  requirements  increase,  as  the  num¬ 
ber  of  I/O’s  multiplies,  and  as  the  clock  frequencies  move 
above  10  GHz,  a  shift  toward  optical  interconnects  probably 
becomes  inevitable.  If  so,  an  additional  thermal  problem  will 
be  realized  because  optical  transmitters  will  need  to  be  inte¬ 
grated  into  the  package  or  perhaps  on-chip.  Research  into 
new  ways  to  improve  the  package  thermal  performance  is 
ongoing,  much  of  which  relies  on  novel  opportunities  offered 
by  nanostructured  materials.  Examples  range  from  high  ther¬ 
mal  conductivity  package  fill  or  epoxy  using  nanotube  com¬ 
posites  to  chip-scale  integrated  thermoelectric  refrigerators. 

The  device  reliability  thermal  problem  refers  most  im¬ 
portantly  to  the  reduction  in  electromigration  lifetime  for 
metallic  interconnect  systems,  in  which  metal  mass  flow  oc¬ 
curs  in  response  to  electrical  current  carried  in  the  intercon¬ 
nects  until  the  metal  line  fails  with  a  high  resistance  condi¬ 
tion.  CMOS  technology  at  dimensional  scales  below  0.25 
jmm  has  significant  performance  limitations  due  to  RC  delays 
in  the  interconnect  network,  so  much  so  that  the  traditional 
Al  based  metallization  has  been  discarded  in  favor  of  Cu  and 
Cu  alloys  and  the  silicon  dioxide  interlayer  dielectric  mate¬ 
rials  are  being  replaced  by  more  exotic  materials  with  lower 
dielectric  permittivity.  In  the  latter  case,  Si02  thin  films  with 
a  relative  permittivity  er^4  is  forecast  to  be  replaced  in  a 
stepwise  fashion  by  materials  with  6r^3  and  soon  after  that 
by  materials  with  er^2.  This  transition  is  being  made  at  a 
rate  limited  by  our  understanding  and  control  of  the  micro¬ 
scopic  and  processing  properties  of  the  films. 

The  advantage  offered  by  these  new  insulator  materials 
has  a  concomitant  set  of  problems,  not  the  least  of  which  is 
an  amplification  of  the  self-heating  problem.  Low  permittiv¬ 
ity  materials  generally  have  low  thermal  conductivity,  the 
consequence  of  which  is  an  increase  in  the  operating  tem¬ 
perature  of  the  metallic  traces,  leading  to  reduced  electromi¬ 
gration  resistance  and  premature  circuit  failure.  Shen102  pub¬ 
lished  a  straightforward  calculation  that  clearly  demonstrates 
this  self-heating  issue.  Examining  a  multilevel  metallization 
idealization  employing  Cu  and  polyimide,  as  a  prototypical 
low  permittivity  material,  the  temperature  increase  in  inter¬ 
connect  system  was  calculated  as  a  function  of  the  number  of 
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FIG.  11.  (Color)  Schematic  representation  of  methyl  doping  of  an  amor¬ 
phous  Si02  network. 

metal  layers  in  the  design.  For  only  a  six  layer  system,  typi¬ 
cal  of  today’s  logic  designs,  the  increase  of  temperature  eas¬ 
ily  exceeded  self-heat  temperature  limits  currently  in  use. 
Future  designs  will  include  more  metal  layers.  Furthermore, 
given  that  logic  circuit  design  employs  a  wide  range  of  local 
metal  densities,  it  is  likely  that  Shen’s  model  underestimates 
the  reliability  jeopardy. 

Given  the  above  considerations,  understanding  the  ther¬ 
mal  conductivity  of  thin  ( ~  <  500  nm)  interlayer  dielectrics 
is  of  considerable  interest.  Using  an  example,  consider  the 
properties  of  carbon  doped  silicon  dioxide  (CDO).103  A  sche¬ 
matic  representation  of  these  materials  is  presented  in  Fig. 
11,  in  which  the  sp 3  network  of  the  amorphous  oxide  is 
opened  up  by  the  introduction  of  chemical  moieties,  such  as 
the  methyl  group  suggested  in  the  figure.  The  thermal  con¬ 
ductivity  of  a  set  of  these  materials  has  been  measured  using 
the  time  domain  thermoreflectance  (TDTR)  method,104  and 
the  results  in  two  cases  independently  validated  using  an 
electrical  resistance  measurement  and  full  scale  finite  ele- 


FIG.  12.  Room  temperature  thermal  conductivities  of  carbon-doped  silicon 
dioxide  (CDO)  films  measured  using  time  domain  thermal  reflectance 
TDTR.  (The  linear  fit  is  merely  a  guide  for  the  eye.) 
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ment  thermal  modeling  method.105  The  results  are  summa¬ 
rized  in  Fig.  12,  in  which  the  measured  thermal  conductivi¬ 
ties  of  thermally  grown  oxide,  plasma  enhanced  TEOS 
oxide,  fluorine  doped  Si02,  a- SiC,  and  various  forms  of 
CDO  are  compared  to  the  measured  oxide  mass  densities, 
determined  using  x-ray  reflectometry.  (The  a  -SiC  samples 
are  actually  a  carbon  rich  CDO  material.)  The  variability  of 
this  measurement  is  estimated  to  be  —5%  determined  by 
testing  multiple  samples  and  over  a  period  of  ^  1  year.  For 
this  family  of  materials  the  introduction  of  carbon  is  the  most 
important  factor  correlated  to  the  drop  in  thermal  conductiv¬ 
ity.  The  preparation  method  and  the  chemical  precursor  pro¬ 
duce  smaller  changes  in  thermal  conductivity.  Carbon  dop¬ 
ing,  in  turn,  reduces  the  resultant  mass  density  of  the  film,  an 
effect  which  is  amplified  by  the  inclusion  of  (CH)X  moieties 
for  x=  1,  2,  and  3  in  the  amorphous  network.106,107  Note  also 
that  F  doping  has  the  same  density  reducing  effect,  albeit 
smaller  in  magnitude. 

The  minimum  thermal  conductivity  theory  of  Cahill  and 
Pohl108  has  been  applied  to  these  results  to  derive  thermal 
conductivities.  This  model  is  an  attractive  option  because  it 
depends  only  on  the  atomic  number  density  and  sound  ve¬ 
locity  of  the  films.  However,  using  the  sound  velocity  de¬ 
rived  from  the  elastic  modulus  measured  separately  using 
nanoindentation  for  each  film,  the  results  obtained  show  only 
a  qualitative  agreement  with  experiment,  with  errors  of 
±  30% .  A  more  successful  approach  would  not  only  close 
the  gap  between  model  and  experiment  but,  for  technological 
relevance,  also  include  the  effects  of  local  chemical  bonding 
and  deposition  process  methodology. 

The  CDO  materials  described  above  are  interesting  in 
that  the  network  is  generally  an  sp3  bonded  one,  in  which 
the  material  density  is  reduced  by  configurational  rearrange¬ 
ment  to  accommodate  the  dopant  materials.  The  resultant 
molecular  cage-like  pores  are  on  the  atomic  scale  and  prob¬ 
ably  do  not  require  explicit  consideration  in  an  elementary 
theoretical  description.  However,  more  aggressive  materials, 
those  providing  a  dielectric  permittivity  near  2,  can  be  pre¬ 
pared  using  the  explicit  introduction  of  nanoscale  porosity 
and  alternative  local  chemistry.  The  manner  in  which  one 
includes  nanoscale  pores  into  a  theoretical  approach  for  ther¬ 
mal  transport  is  not  yet  precise.  For  example,  the  description 
of  the  transition  from  “cage-like”  to  nanoporosity  needs  to 
be  clarified.  Perhaps  more  importantly,  the  pore  sizes  of  in¬ 
terest,  which  are  a  few  nanometers  so  as  not  to  limit  the 
mechanical  strength  of  the  films,  are  comparable  to  the  rel¬ 
evant  thermal  phonon  wavelengths,  which  renders  con¬ 
tinuum  approximations  suspect.  The  effect  of  porosity  on  the 
room  temperature  thermal  conductivity  of  xerogel,  a  nanopo- 
rous  form  of  a- Si02,  was  recently  studied  by  Hu  and 
co-workers;109  they  proposed  phenomenological  two-phase 
mixture  models  that  usefully  describe  the  dependence  of  the 
conductivity  on  porosity. 

To  explore  the  effect  of  porosity  on  the  thermal  transport 
in  lower  permittivity  amorphous  materials,  Costescu  et  al.no 
have  investigated  the  thermal  conductivity  of  hydrogen- 
silsesquioxane  (HSQ)  films  in  the  temperature  range  80<  T 
<400  K  using  the  3co  method.13,111  HSQ  is  a  spin-on  mate¬ 
rial  with  dielectric  permittivities  near  to  2.  Figure  13  shows 
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FIG.  13.  Comparison  of  data  for  hydrogen-silsesquioxane  HSQ  and  Si02 
(filled  circles)  to  calculations  of  the  minimum  thermal  conductivity  (dashed 
lines)  for  each  material.  The  mass  density  of  the  “flowable  oxide”  sample, 
FOx-A,  is  p=  1.5  gem-3;  for  the  “extra-low-/:”  sample,  XLK-6,  p 
=  0.90  g  cm-3. 

the  results  of  thermal  conductivity  measurements  on  a  vari¬ 
ety  of  HSQ  samples  that  differ  primarily  by  processing  con¬ 
ditions,  which  is  expected  to  modulate  the  porosity  of  result¬ 
ant  films.  Remarkably,  the  temperature  dependence  of  the 
data  is  nearly  identical  for  all  samples,  i.e.,  the  conductivity 
of  HSQ  is  suppressed  relative  to  Si02  by  a  constant  factor 
that  is  independent  of  temperature  in  the  temperature  range 
of  the  experiments,  80<r<400  K,  and  a  significant  differ¬ 
ence  in  porosity  does  not  appear  to  modify  the  character  of 
thermal  transport. 

Thus,  the  nearly  identical  temperature  dependence  of  the 
thermal  conductivities  of  Si02 ,  FO* ,  and  XLK  suggests  that 
the  spectrum  of  heat-carrying  vibrational  modes  is  nearly  the 
same  in  all  three  materials.  It  is  thus  attractive  to  conclude 
that  heat  transport  in  these  materials  is  intrinsically  local  and 
that  porosity  does  not  play  a  role  except  to  reduce  average 
density.  However,  this  is  the  basic  physical  picture  of  the 
minimum  thermal  conductivity  theory  applied  above  for  the 
doped  oxides.108  The  results  of  this  theory  applied  to  the 
HSQ  materials  are  presented  in  Fig.  13,  where  it  is  evident 
that  with  increasing  porosity,  there  is  an  increasingly  poor 
match  between  the  data  and  the  model.  Effective  medium 
theory  (Refs.  112  and  113)  in  which  the  porous  solid  is 
treated  as  a  composite  material  of  a  matrix  punctuated  by 
voids  offers  a  path  to  better  understand  these  materials.110 

B.  Multilayers  and  superlattices 

Multilayers  and  superlattices  are  films  that  contain  alter¬ 
nating  layers  of  two  different  materials  stacked  upon  each 
other.  In  multilayers  the  films  may  be  amorphous  or  poly¬ 
crystalline  while  in  superlattices  the  films  are  single  crystal. 
In  this  section,  we  concentrate  on  the  thermal  conductivity  in 
the  direction  perpendicular  to  the  plane  of  the  film  because 


this  is  the  property  accessible  to  the  most  commonly  used 
experimental  techniques:  time-domain  thermoreflectance,  see 
Sec.  V  A,  and  the  3co  method  for  thin  films.111,13  (When  the 
width  of  the  heater/thermometer  line  in  the  3co  method  is 
comparable  to  the  film  thickness,  the  3a>  method  is  also  sen¬ 
sitive  to  the  thermal  conductivity  in  the  plane  of  the  film.114 
Data  collected  using  combinations  of  narrow  and  wide  metal 
lines  can  be  used  to  separate  the  in-plane  and  cross-plane 
conductivities.114) 

Phonon  transport  in  superlattices  and  multilayers  can  be 
divided  into  two  categories,  namely:  (i)  €</?;  and  (ii)  € 
>b ,  where  €  is  the  phonon  mean  free  path  and  b  is  the 
superlattice  period.  Consider  the  first  regime.  There  have 
been  a  number  of  studies  of  thermal  transport  in  multilayer 
films12,14,115  with  heat  transport  perpendicular  to  the  plane  of 
the  film.  The  general  conclusion  from  these  studies  is  that  if 
the  thermal  conductivity  of  a  thin  layer  of  disordered  mate¬ 
rial  is  independent  of  layer  thickness  b ,  i.e.,  if  phonons  with 
mean-free-paths  comparable  to  b  make  a  negligible  contri¬ 
bution  to  heat  transport,  and  that  the  micro  structure  of  the 
thin  layer  is  approximately  independent  of  b,  the  measured 
resistance  of  a  multilayer  sample  is  the  series  sum  of  the 
thermal  resistance  of  the  individual  layers  and  the  thermal 
resistance  of  the  interfaces.13,116  Thus,  a  simple  way  of  pre¬ 
dicting  the  thermal  conductivity  is  to  add  up  the  bulk  plus 
interfacial  resistances  of  all  the  layers. 

Figure  1  shows  the  data  and  predictions  for  the  thermal 
conductance  of  a  single  interface.  The  values  range  from  3  to 
200  MW  m-2  K- 1  at  room  temperature  while  the  predictions 
of  the  diffuse  mismatch  model  (DMM)  tend  to  be  on  the 
order  of  500  MWnrG2  K-1.  While  this  discrepancy  contin¬ 
ues  to  exist,  in  many  cases  the  interfacial  resistance  is  neg¬ 
ligible  compared  to  the  bulk  resistance  of  a  single  layer.  This 
often  occurs  when  the  films  are  amorphous.  In  such  cases, 
the  thermal  conductivity  of  multilayer  films  is  found  to  be 
independent  of  period  b .  Studies  of  multilayers  of  disordered 
oxides  have  failed  to  reveal  a  significant  effect  of  interfaces: 
the  thermal  conductivity  of  Zr02 :  Y203  /Si02  (Ref.  117)  and 
Zr02/Y203  (Ref.  12)  multilayers  are  almost  independent  of 
layer  thickness  with  layer  thickness  as  small  as  4  nm.  Appar¬ 
ently,  the  lattice  vibrations  of  these  oxides  are  sufficiently 
alike  that  the  thermal  conductance  is  too  large  to  be  observed 
by  this  approach.  Multilayer  coatings115  do  not,  therefore, 
appear  to  be  a  promising  route  for  lowering  the  conductivity 
of  thermal  barriers.  Significant  decreases  in  conductivity 
have  been  observed  recently  in  nanocrystalline  Zr02:Y203 
(Ref.  118)  and  attributed  to  the  finite  thermal  conductance  of 
grain  boundaries. 

Theoretical  and  experimental  studies  on  thermal  trans¬ 
port  in  epitaxial  superlattices  have  been  conducted  for  more 
than  two  decades.  In  most  cases,  phonon  transport  in  super¬ 
lattices  fall  into  the  second  regime,  i.e.,  (>b.  Hence,  both 
wave  and  particle-like  phonon  effects  are  observed,  offering 
not  only  richness  in  physics  but  also  difficulties  in  predic¬ 
tions.  The  origins  of  this  work  can  be  traced  back  to  Naray- 
anamurti  et  al.'}19  they  studied  monochromatic  phonon  trans¬ 
port  across  lattice  matched  GaAs/AlGaAs  superlattices  and 
showed  the  existence  of  a  phonon  band  gap.  This  is  a 
coupled  behavior  resulting  from  interference  of  waves  re- 
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TABLE  I.  Experimental  research  on  thermal  transport  in  semiconductor  superlattices.  Unless  otherwise  noted,  the  comments  refer  to  cross-plane  transport. 
AM  stands  for  “acoustic  mismatch.” 


Lowest  thermal 


Reference 

Materials 

Effects  considered 

Temperature  behavior 

Period  behavior 

conductivity 

124 

GaAs/AlAs 

not  applicable 

not  applicable 

in-plane:  max  at  200  A 

in-plane:  12  Wm-1  K_1 

125 

GaAs/AlAs 

interface  roughness 

in-plane:  decrease  190-450  K 

(700  A/700  A  only) 

in-plane:  ^40 

133 

GaAs/AlAs 

AM,  band  gaps,  roughness 

decrease  100-400  K 

increase  10-400  A 

4 

122 

Si/Ge 

dislocations 

increase  to  150  K 

maximum  at  70  A 

«1.6 

67 

Si/Sio.7iGe0.29 

AM 

not  applicable 

50  A/50  A  only) 

22,  in-plane:  54 

123 

Si/Ge 

doping,  period, 
dislocations 

increase  80-150  K, 
constant  to  330  K 

decrease  b  =  40-  90  A 

slight  decrease  to  b=  140  A 

1.2 

127 

Si/Sio  yGeg  3 

AM 

increase  100-300  K 

increase  for  b  =  45-300  A 
in-plane:  max  shifts  with  period 

10 

134 

Sio.84Geo.i6  /  Sio.76Geo.24 

AM 

increase  50-180  K, 
constant  to  300  K 

constant  for  b  =  61,  100,  133  A 

8.4 

128 

InAs/AlSb 

growth  temp,  annealing 

maximum  at  150  K 

not  applicable 

2.5 

135 

InP/InGaAs 

thickness  ratio 

increase  50-100  K 

not  applicable 

7 

132 

Bi2Te3  /Sb2Te3 

doping,  tunneling, 
localization 

not  applicable 

lattice:  minimum  at  40  A 

0.2 

12 

PbTe/PbSe 

not  applicable 

decrease  200-400  K 

40  A  only 

1.6 

fleeted  from  multiple  interfaces.  For  this  to  occur,  the  mean 
free  path  of  phonons  must  span  multiple  interfaces  (>  10). 
In  this  regime,  the  phonon  dispersion  relation  is  modified 
and  zone  folding  occurs  resulting  in  multiple  phonon  band 
gaps.120  Several  effects  can  result  from  this  modification  of 
the  phonon  dispersion.  First,  the  group  velocities  of  phonons 
are  reduced  significantly,  especially  for  higher  energy  acous¬ 
tic  phonons.  Second,  because  co-k  relations  are  modified, 
there  are  many  more  possibilities  for  conservation  of  crystal 
momentum  and  energy  involved  in  normal  and  umklapp 
scattering.121  Hence,  the  scattering  rate  is  increased. 

In  the  regime  when  the  mean  free  path  is  not  large 
enough  for  phonon  band  gaps  to  occur,  single  interface  ef¬ 
fects  can  also  play  a  significant  role.  For  example,  if  the  two 
materials  in  the  superlattice  have  large  mismatch  in  the  pho¬ 
non  dispersion  relations,  phonons  in  certain  frequency  range 


cannot  propagate  to  the  neighboring  layer  unless  there  are 
mode  conversions  at  the  interface.  In  addition  to  these  ef¬ 
fects,  interfaces  between  two  different  materials  with  differ¬ 
ent  lattice  constants  can  contain  dislocations  and  defects. 
These  can  also  scatter  phonons  and  reduce  thermal  conduc¬ 
tivity.  Finally,  depending  on  the  type  of  processing,  there 
could  be  both  physical  roughness  as  well  as  alloying  at  the 
interface,  which  can  also  influence  phonon  transport.  Hence, 
it  is  clear  that  there  are  many  different  effects  that  contribute 
to  the  resistance  to  heat  flow. 

To  dissect  the  contributions  of  these  effects,  several  ex¬ 
perimental  and  theoretical  studies  have  been  performed. 
Table  I  lists  the  experimental  work  done  in  this  area  while 
Table  II  lists  the  theoretical  and  computational  ones.  Most 
experimental  studies  have  used  superlattices  of  Si,  Ge,  and 
their  alloys  as  well  as  those  of  III-V.  Table  I  also  shows  the 


TABLE  II.  Theoretical  research  on  thermal  transport  in  semiconductor  superlattices.  AM  stands  for  “acoustic  mismatch.” 


Reference 

Superlattice 

Model 

Effects  considered 

Comparison  with  experiment 

121 

general 

phenomenological 

minibands 

conductivity  increases  with  period 

136 

Si/Ge 

sc  lattice  dynamics 

AM,  phonon  confinement,  minibands 

room  temperature  reduction 
agrees  well  with  Si/Ge 

137-139 

Si/Ge; 

GaAs/AlAs 

Boltzmann  transport 
equation 

AM,  inelastic  scattering,  roughness 

see  references 

140 

Si/Ge 

fee  lattice  dynamics 

minibands,  velocity 

trend  does  not  agree  with  experimental  results 

141 

general 

transfer  matrix 

phonon  interference,  tunneling 

poor  agreement  with  experiment 

142 

GaAs/AlAs 

lattice  dynamics 

velocity  reduction 

poor  agreement  with  experiment 

143 

Si/Ge 

valence-force  potential 

velocity  reduction,  spectra  mismatch 

periods  too  thin  for  comparison 

144 

Si/Ge  &  GaAs/AlAs 

fee  lattice  dynamics 

minibands,  tunneling 

good  agreement  for  cross-plane; 
poor  agreement  in-plane 

120 

general 

wave  theory 

mini-bands,  interface  scattering 

trend  agrees  with  lattice  conductivity  of  Bi2Te3  /Sb2Te3 

145 

Si/Ge 

molecular  dynamics 

all  classical  effects 

good  agreement  (two  points)  with  experiment 

146 

general 

phenomenological 

particle-wave  crossover 

high  TC  at  ultra-short  periods;  model 
cannot  explain  trend 

147 

SiGe  quantum  dot 

continuum 

quantum  dot  scattering 

not  applicable 

132 

Bi2Te3  /Sb2Te3 

Boltzmann  transport 
equation 

AM,  interface  roughness, 
localization 

trend  does  not  agree  with  data 

148 

Lennard-Jones  solids 

molecular  dynamics 

all  classical  effects,  strain 

not  applicable 

149 

fee 

molecular  dynamics 

all  classical  effects 

not  applicable 
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FIG.  14.  Thermal  conductivity  of  Si/Ge  (top)  (Ref.  122)  and  Si/Si0.7Ge0.3 
(bottom)  (Ref.  127)  superlattices;  the  curves  are  labeled  by  the  superlattice 
period  in  A.  The  thermal  conductivity  of  Si/Ge  superlattices  (top)  increases 
with  period  thickness  in  the  region  of  30-60  A,  but  then  drops  dramatically 
for  period  higher  than  140  A,  presumably  due  to  presence  of  misfit  disloca¬ 
tions  and  defects  beyond  a  critical  thickness.  The  thermal  conductivity  of 
Si/Si0  7Ge0  3  with  thickness  ratio  of  2:1  showed  that  as  the  period  was  de¬ 
creased  from  300  A  (100  periods  in  3  jul m  thick  him)  to  45  A  (667  periods 
in  3  jmm  thick  him),  the  thermal  conductivity  decreased  steadily  and  ap¬ 
proached  the  effective  alloy  which  is  Si0  9Ge0  x . 


effects  considered  (interfacial  roughness,  phonon  band  gaps, 
dispersion  mismatch,  etc.)  and  the  observed  dependence  of 
thermal  conductivity  on  temperature  and  superlattice  period. 
The  data  on  Si  and  Ge  superlattices  suggest  that  for  very 
small  period,  when  the  layer  thickness  is  less  than  the  critical 
one  for  the  extension  of  misfit  disclocations,  the  thermal  con¬ 
ductivity  of  superlattices  tend  to  decrease  with  decreasing 
layer  thickness.122  As  the  layer  thickness  is  increased  beyond 
10  nm,  the  thermal  conductivity  drops,  presumably  due  to 
disorder  created  by  extensive  plastic  deformation  [see  Fig. 
14(a)].  This  is  similar  to  the  observations  of  BorcaTascuic 
et  al. 123 

Interface  disorder  introduces  diffuse  scattering  but  inter¬ 
face  disorder  in  a  superlattice  is  typically  modest.  It  is  lim¬ 
ited  to  interface  roughness  (a  finite  density  of  interface  steps) 
and  substitutional  alloying  at  the  interfaces  caused  by  surface 
segregation  during  growth.  Unfortunately,  this  physical  and 
chemical  interface  roughness  in  an  epitaxial  superlattice  is 
sensitive  to  the  material,  growth  method,  growth  tempera¬ 
ture,  and  deposition  rate — and  is  notoriously  difficult  to  char¬ 
acterize.  Furthermore,  asymmetry  of  the  interface  is  often 
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FIG.  15.  Selected  data  for  the  through-thickness  thermal  conductivity  of 
superlattices  with  bilayer  periods  of  «5nm;  data  for  GaAs-AIAs  [open 
circles,  5.67  nm  period  (Ref.  126)]  were  measured  by  picosecond  thermore¬ 
flectance;  data  for  Si-Ge  [open  triangles,  5  nm  period  (Ref.  122);  filled 
circles,  4.4  nm  period  (Ref.  123)]  and  InAs-AISb  [open  diamonds,  6.5  nm 
period  (Ref.  128)  were  measured  using  the  3 co  method.  Data  for  Si0  85Ge0 15 
(Ref.  122)  and  Gao4Al06As  (Ref.  129)  alloys  are  included  for  comparison. 


pronounced,  i.e.,  the  growth  of  component  A  on  the  surface 
of  component  B  creates  a  different  interface  than  the  growth 
of  B  on  A.  Strained-layer  superlattices  such  as  Si-Ge  can 
harbor  high  densities  of  crystalline  defects  when  the  layer 
thickness  exceeds  the  critical  thickness  for  the  extension  of 
misfit  dislocations;  growth  on  relaxed  buffer  layers  reduces 
the  density  of  threading  dislocations  but  the  misfit  density 
will  still  be  large  when  the  critical  thickness  is  exceeded.  The 
agreement  between  the  data  of  Cahill  et  al}22  and  BorcaTas¬ 
cuic  et  al 123  seem  to  support  this  idea.  The  interesting  ob¬ 
servation  that  they  both  made  was  that  the  superlattice  ther¬ 
mal  conductivity  was  lower  than  that  of  the  alloy  with  the 
same  mass  ratios  (see  Fig.  15).  Hence,  interface  scattering 
alone  is  probably  insufficient  to  explain  this  observation. 
Perhaps  this  may  be  due  to  zone  folding  which  reduces  pho¬ 
non  group  velocities  and  increases  scattering  rates  and 
thereby  reduces  thermal  conductivity.  Thermal  transport  in 
GaAs  and  AlAs  superlattices  have  also  been  studied 
widely.66,124-126  These  studies  have  shown  that  as  the  inter¬ 
face  density  increases,  the  thermal  conductivity  decreases, 
clearly  showing  that  interface  resistance  is  important. 

More  recently,  Huxtable  et  al.121  have  shown  that  for 
superlattices  of  Si  and  Si0  7Ge0  3 ,  where  the  critical  thickness 
is  much  larger  than  that  for  Si  and  Ge  superlattices,  the 
thermal  conductivity  scaled  almost  linearly  with  interface 
density  supporting  the  data  of  Cahill  et  al 122  and 
BorcaTascuic123  [see  Fig.  14(b)].  As  the  interface  density  in¬ 
creased,  the  thermal  conductivity  approached  that  of  the  al¬ 
loy,  but  was  never  lower  than  that  of  the  alloy,  which  is  in 
contrast  with  previous  observations.  Clearly,  the  mismatch  in 
acoustic  impedance  between  Si  and  Ge  is  more  than  that  of 
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Si  and  Si0  7Ge0  3 .  Hence,  it  is  possible  that  thermal  conduc¬ 
tivities  below  alloy  values  can  be  found  only  when  there  is 
large  mismatch  in  acoustic  impedance  that  creates  significant 
band  gaps  in  the  phonon  dispersion  relation. 

There  have  been  two  reports  of  thermal  conductivity  of 
Bi2Te3  and  Sb2Te3  superlattices.  These  are  of  interest  in  ther¬ 
moelectric  applications  where  the  reduction  of  thermal 
conductivity  has  been  a  major  challenge.130 
Venkatasubramanian131  showed  as  the  period  thickness  in¬ 
creased,  the  room  temperature  thermal  conductivity  reached 
a  minimum  for  a  period  thickness  5  nm,  and  this  value  was 
lower  than  the  thermal  conductivity  of  the  corresponding  al¬ 
loy.  It  was  suggested  that  this  was  due  to  phonon  band-gap 
effects.  However,  Touzelbaev  et  al.132  did  not  find  any  sig¬ 
nificant  change  in  thermal  conductivity  when  the  superlattice 
period  was  varied.  Hence,  it  remains  inconclusive  as  to 
whether  phonon  band-gap  effects  are  indeed  present  in  this 
system. 

Table  II  lists  the  theoretical  and  computational  work  on 
thermal  transport  in  superlattices.  Three  approaches  are  gen¬ 
erally  adopted.  The  first  is  based  on  wave  propagation,  the 
second  on  particle  transport  that  can  be  modeled  by  Boltz¬ 
mann  transport  theory,  and  the  third  relies  on  molecular  dy¬ 
namics. 

From  the  experimental  and  theoretical  studies,  it  is  clear 
that  the  mechanism  of  phonon  transport  in  epitaxial  superlat¬ 
tices  is  not  well  understood  and  it  is  difficult  to  predict  the 
value  of  thermal  conductivity.  However,  a  few  trends  do 
seem  to  appear.  First,  if  interface  resistance  is  to  play  a  major 
role,  the  acoustic  mismatch  between  the  two  films  of  the 
superlattice  must  be  sufficiently  large,  >1.1.  For  incoherent 
phonons,  acoustic  mismatch  gives  rise  to  interface  reflection 
while  for  coherent  phonons,  phonon  band  gaps  are  formed 
that  reduce  the  group  velocities.  Second,  interface  roughness 
and  alloying  can  also  scatter  phonons  and  thereby  increase 
interfacial  resistance.  Third,  in  some  cases,  the  phonon  ther¬ 
mal  conductivity  can  be  reduced  below  the  alloy  limit,  which 
may  occur  due  to  band-gap  formation.  Although  calculations 
predict  this,  experiments  have  not  provided  unmistakable 
proof.  The  image  that  seems  to  emerge  from  all  these  studies 
is  the  following.  To  reach  the  minimum  thermal  conductivity, 
one  perhaps  needs  to  take  different  approaches  to  block  pho¬ 
non  transport  in  different  parts  of  the  phonon  spectrum.  For 
example,  high-frequency  phonons  can  be  blocked  by  alloy 
scattering  since  the  wavelengths  are  generally  on  the  order  of 
a  few  atomic  spacings.  Low  frequency  phonons  are  generally 
unaffected  by  alloy  scattering  since  the  scattering  is  gener¬ 
ally  in  the  Rayleigh  regime.  However,  they  may  dominate 
energy  transport  in  an  alloy.  For  such  phonons,  superlattice 
formation  would  lead  to  miniband  formation  since  their 
mean  free  paths  are  generally  much  larger  than  for  high- 
frequency  phonons.  Hence,  what  would  indeed  be  interesting 
is  to  study  thermal  transport  in  superlattices  of  alloys  that 
have  sufficient  acoustic  impedance  mismatch. 

C.  Superlattices:  Theory  and  computation 

Thermal  transport  in  semiconductor  superlattices  can  be 
either  parallel  or  perpendicular  to  the  plane  of  the  layers.  The 
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initial  measurements  and  theories  were  for  transport  parallel 
to  the  layers.  The  simple  viewpoint  is  that  heat  transport 
should  be  quite  efficient  parallel  to  the  layers.  If  the  inter¬ 
faces  are  atomically  smooth,  then  phonons  would  specularly 
reflect  from  these  interfaces.  Each  layer  becomes  a  phonon 
waveguide  which  efficiently  channels  heat  along  each  layer. 
The  initial  measurements  of  thermal  conductivity  parallel  to 
the  layers  found  values  smaller  than  bulk  values  by  a  factor 
of  four.  Yao124  found  in  GaAs/AlAs  superlattices  that  the 
thermal  conductivity  was  similar  to  that  of  a  random  alloy  of 
the  two  materials.  This  was  explained  in  Ref.  150  as  due  to 
the  interface  boundaries  not  being  perfectly  planar.  Atomic 
scale  defects  at  the  interface  are  efficient  scattering  centers 
for  the  phonons. 

There  is  much  more  interest  in  thermal  transport  perpen¬ 
dicular  to  the  plane  of  the  layers,  which  is  called  the  “cross¬ 
plane”  direction.  In  some  cases  experiments  show  a  thermal 
conductivity  significantly  lower  than  that  of  the  random  alloy 
of  the  two  materials  in  the  superlattice.122,131,133  The  simple 
explanation  for  this  reduction  is  that  the  superlattice  acts  as  a 
Fabrey-Perot  interferometer,  and  allows  only  a  few  phonon 
frequencies  to  pass  through.  This  feature  was  demonstrated 
by  phonon  transmission  spectroscopy.119 

There  have  been  many  calculations  of  thermal  conduc¬ 
tivity  in  the  cross  plane  direction.  They  can  be  divided  into 
three  groups.  The  largest  number  of  calculations  are  based  on 
the  standard  formula  for  the  thermal  conductivity 

^  f  d3q  _  I  dnB[a>x(q)]\ 

K  =  h ?  J  (2^VM  ^(qK(q)( - df - )* 

(8) 

where  uzX(q)  is  the  velocity  of  the  phonon,  rx(q)  is  the 
lifetime,  and  nB[o)x( q)]  is  the  Boson  occupation  function. 
The  above  formula  is  derived  from  the  Boltzmann  transport 
equation  (BTE).  Calculation  by  Kato  et  al,151  Hyldgaard  and 
Mahan,136  Tamura  et  al.,140  Simkin  and  Mahan,120  Bies 
et  al.}42  and  Kiselev  et  al.}43  all  just  calculate  the  actual 
phonon  modes  of  the  superlattice  and  evaluate  the  above 
formula.  Usually  the  lifetime  rx(q),  or  else  the  mean-free- 
path  €Xq  is  selected  to  be  the  same  as  in  the  homogeneous 
material.  These  calculations  show  the  thermal  conductivity  in 
the  superlattice  is  reduced  significantly,  sometimes  by  a  fac¬ 
tor  of  ten,  compared  to  that  of  the  constituents  of  the  super¬ 
lattice.  This  large  reduction  agrees  with  the  experimental 
findings. 

Simkin  and  Mahan120  used  a  complex  wave  vector  in 
computing  the  superlattice  modes.  The  imaginary  part  of  the 
wave  vector  was  a  phenomenological  way  of  including  the 
mean-free-path  into  the  calculation  of  the  phonon  modes. 
This  calculation  showed  a  minimum  in  the  thermal  conduc¬ 
tivity  as  a  function  of  superlattice  period. 

There  have  been  several  calculations  of  the  thermal  con¬ 
ductivity  of  the  superlattice  which  have  solved  the  Boltz¬ 
mann  transport  equation.  That  of  Ren  and  Dow121  did  not 
predict  a  large  reduction  of  the  thermal  conductivity,  and 
disagrees  with  experiments.  Several  calculations  from 
Chen’s  group138,152  agree  much  better  with  experiment. 

Recently,  Daly  and  co-workers149,153  used  large-scale 
(1.3X106  atoms)  molecular  dynamics  simulations  to  calcu- 
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late  the  thermal  conductivity  of  superlattices. 

There  have  been  very  few  calculations  of  thermal  con¬ 
ductivity  in  the  cross  plane  direction  which  included  heat 
transfer  by  electrons  as  well  as  phonons.  Bartkowiak  and 
Mahan’s  modeling146  was  semiclassical  and  assumed  that 
both  the  electrons  Te(x)  and  phonons  Tp(x)  had  a  local  tem¬ 
perature.  The  electrons  and  phonons  can  exchange  energy  in 
the  bulk,  and  at  the  interfaces.  They  found  a  wide  variety  of 
temperature  distributions  depending  upon  the  various  bound¬ 
ary  resistances  for  electrons  and  phonons.  None  of  this  be¬ 
havior  has  been  verified  or  refuted  by  experimental  measure¬ 
ments.  All  of  their  calculations  depend  on  the  concept  of  a 
local  temperature. 


V.  NEEDS  AND  LIMITS  OF  METROLOGY 
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The  goal  of  successful  commercial  implementation  of 
nanoscale  materials  places  a  burden  on  the  metrology  of 
these  materials  and  devices  that  extend  beyond  the  charac¬ 
terization  of  the  fundamental  materials  properties.  This  bur¬ 
den  involves  not  only  the  small  physical  dimension  of  the 
material  but  also  the  general  accuracy  of  the  measurement. 
To  illustrate  this,  perhaps  the  most  familiar  examples  derive 
from  the  Si  semiconductor  industry.  Here,  thin  film  metrol¬ 
ogy  is  becoming  increasingly  challenging  as  Si  semiconduc¬ 
tor  device  technology  scales  toward  smaller  dimensions  and 
the  corresponding  manufacturing  process  complexity  grows. 
Traditional  thin  film  properties  such  as  thickness  and  intrin¬ 
sic  stress  are  being  augmented  by  an  increasing  number  of 
other  physical  parameters,  which  together  must  be  consid¬ 
ered  in  the  development  and  manufacturing  monitoring  plan. 
Such  parameters  include  density,  porosity,  thermal  conduc¬ 
tivity,  grain  texture,  surface  chemical  residue,  and  bulk  and 
interface  cohesion.  The  introduction  of  more  complex  de¬ 
vices  using  new  materials  and  process  methods  (e.g.,  Si  al¬ 
loys,  polymers,  and  3D  structures)  also  contribute  to  the 
measurement  challenge.  Effective  control  of  these  properties 
often  becomes  difficult  as  the  device  scaling  continues.  For 
example,  the  control  of  the  physical  thickness  and  composi¬ 
tion  profile  of  the  MOS  gate  dielectric  is  generally  of  critical 
importance  to  the  resultant  device  parameters.  But  with  this 
film  currently  having  a  required  thickness  of  less  than  1  nm, 
the  challenge  is  to  control  thickness  and  composition  varia¬ 
tion  on  the  atomic  scale. 

Another  factor  to  consider  is  that  the  ability  to  measure 
many  properties  is  becoming  less  and  less  feasible  at  the 
nanoscale  limit.  A  well  known  example  is  the  measurement 
of  dopant  distributions  under  a  gate  as  the  gate  widths  move 
below  the  50  nm  scale  and  approach  the  “sparse  dopant” 
limit.154  As  a  result,  an  increasing  number  of  traditional  mea¬ 
surement  methods  are  becoming  insensitive  or  unsuitable  for 
use  due  to  the  limitations  of  the  physics  employed.  To  ad¬ 
dress  these  difficulties,  not  only  is  new  instrumentation  being 
sought,  but  also  industry  is  relying  increasingly  on  materials, 
device,  and  thermal-mechanical  modeling  to  infer  properties 
that  cannot  be  directly  measured.  This  approach,  in  turn,  re¬ 
quires  that  the  fundamental  properties  be  precisely  and  accu¬ 
rately  measured  on  test  systems,  designed  to  enable  such 
measurements. 


FIG.  16.  Schematic  of  the  TDTR  and  picosecond  acoustics  apparatus  in  use 
at  the  Laser  Facility  of  the  Frederick  Seitz  Materials  Research  Laboratory. 

The  consequence  of  the  above  considerations  for  the 
problem  of  heat  transport  in  nanostructures  has  been  the  in¬ 
troduction  of  a  wide  variety  of  novel  instrumentation  well 
suited  for  measurement  of  thin  films  and  structures.  Several 
promising  approaches  have  been  introduced  recently  or  are 
currently  under  development.  Three  examples  are  described 
below. 

A.  Time-domain  thermoreflectance 

The  first  thermal  transport  experiments  using  picosecond 
thermoreflectance,  alternatively  referred  to  as  time-domain 
thermoreflectance  (TDTR),  were  reported  in  Refs.  155  and 
156.  TDTR  provides  a  direct  method  for  measuring  heat  dif¬ 
fusion  on  nanometer  length  scales.  Most  engineering  materi¬ 
als  have  thermal  diffusivities  D  in  the  range  0.005 <D 
<lcm2s_1  and  at  f  =  200ps,  heat  diffusion  lengths  l 
=  yjDt  are  in  the  range  10</<  140  nm.  Therefore,  picosec¬ 
ond  time  resolution  offers  nanometer- scale  depth  resolution 
and,  more  specifically,  picosecond  thermoreflectance  can  iso¬ 
late  the  effects  of  interface  conductance  from  the  thermal 
conductivity  of  a  thin  layer.  By  contrast,  thermoreflectance 
measurements  using  longer  pulse  length  Q -switched  lasers157 
and  the  3co  method13  cannot  distinguish  between  the  thermal 
conductivity  of  a  film  and  the  thermal  conductance  of  its 
interfaces. 

An  example  of  TDTR  experimental  apparatus  is  pre¬ 
sented  in  Fig.  16.  Picosecond  time  resolution  is  obtained 
using  mode-locked  lasers  that  produce  a  series  of  <  1  ps 
pulses  at  a  repetition  rate  of  ^76  MHz.  Extremely  high  time 
resolution  is  achieved  by  splitting  the  laser  output  into  two 
beam  paths,  a  “pump”  beam  and  a  “probe”  beam,  and  ad¬ 
justing  the  relative  optical  path  lengths  with  a  mechanical 
delay  stage.  If  the  sample  under  investigation  is  metallic,  or 
can  be  coated  with  a  thin  metal  film,  then  a  small  fraction  of 
the  energy  from  each  pulse  in  the  pump  beam  produces  a 
sudden  jump  in  temperature  of  ~  3  K  near  the  surface  of  the 
sample.  The  decay  of  the  near-surface  temperature  is  then 
interrogated  by  the  reflected  energy  of  the  series  of  pulses  in 
the  probe  beam.  Determination  of  the  thermal  conductivity 
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and  interfacial  thermal  resistance  is  made  by  comparing  the 
experimental  cooling  curve  to  the  theoretical  model  and  op¬ 
timizing  the  free  parameters.  For  all  but  the  shortest  times, 
the  temperature  is  homogeneous  through  the  optical  absorp¬ 
tion  depth  of  the  thin  metal  film  and  A/?!  (7) 
=  (dR/dT)ATl(t).  The  surface  temperature  AT^f)  can  be 
accurately  calculated  using  one-dimensional  heat  flow  since 
the  thermal  diffusion  length  is  usually  much  smaller  than  the 
radius  of  the  focused  pump  beam.133 

In  a  modulated  pump -probe  experiment,  the  differences 
in  reflected  probe  intensity  caused  by  the  pump  pulse  appear 
at  the  modulation  frequency  of  the  pump  beam  and  are  ex¬ 
tracted  with  lock-in  detection  but  interpretation  of  this 
lock-in  signal  is  not  straightforward  as  the  time  delay  be¬ 
comes  a  significant  fraction  of  the  separation  between  pulses: 
Capinski  and  co-workers133  discuss  calculations  of  the  in- 
phase  and  out-of-phase  signals  of  the  lock-in  amplifier  given 
the  time-domain  thermal  response  of  the  sample;  Cahill 
et  al. 1  provide  a  calculation  starting  from  the  frequency- 
domain  response. 

Despite  the  promise  of  these  methods,  few  quantitative 
and  systematic  studies  of  thermal  properties  of  materials 
have  used  picosecond  optical  techniques.10,11,133,156,158,159 
Experiments  by  Maris  and  co-workers  on  interface  thermal 
conductance10  and  superlattice  conductivity133  are  the  most 
complete.  Data  for  the  thermal  decay  at  short  times  t 
<  50  ps  are  difficult  to  interpret  quantitatively  because  (i) 
hot-electrons  can  deposit  energy  beyond  the  optical  absorp¬ 
tion  depth,  (ii)  the  temperature  dependence  of  the  complex 
index  of  refraction  dn/dT  is  typically  unknown,  (iii)  the  dif¬ 
fusion  equation  fails  on  length  scales  comparable  to  the 
mean-free-paths  of  the  dominant  energy  carriers,  and  (iv)  the 
assumption  of  equilibrium  between  phonons  and  electrons  is 
not  always  valid.  At  long  times  £>500  ps,  measurements  are 
often  plagued  by  shifts  in  the  focal  plane  and  beam  overlap 
created  by  large  displacements  of  the  delay  stage.  (Capinski 
and  Maris104  describe  an  optical  design  that  improves  accu¬ 
racy  at  long  delay  times;  the  out-of-phase  signal  of  the 
lock-in  amplifier  can  also  be  used  to  correct  systematic 
errors.1)  It  is  likely,  however,  that  improved  performance  and 
dropping  costs  of  solid-state  mode-locked  lasers  will  expand 
applications  of  picosecond  thermoreflectance.  Reliable,  low- 
noise  Ti:  sapphire  lasers  that  produce  100  fs  pulses  are  com¬ 
mercially  available  and  can  be  operated  with  only  a  modest 
knowledge  of  the  technology  of  ultrafast  lasers. 

Over  the  longer  term,  improved  laser  performance  will 
also  make  more  sophisticated  experiments  possible.  For  ex¬ 
ample,  the  pump  and  probe  light  pulses  can  be  applied  to 
localized  regions  of  the  sample  through  the  use  of  near  field 
optics.  With  this  technique,  it  should  be  possible  to  study  the 
lateral  flow  of  heat  in  a  thin  film,  instead  of  just  the  flow  of 
heat  through  the  film.  In  addition,  the  accuracy  of  TDTR 
measurements  will  be  enormously  enhanced  if  compact  laser 
sources  that  can  be  tuned  over  a  wide  wavelength  range  be¬ 
come  available.  The  probe  wavelength  could  then  be  chosen 
to  match  a  peak  in  the  thermoreflectance  of  the  metal  film. 
This  would  increase  signal  to  noise  in  the  measurement  pro¬ 
cess,  reduce  the  time  needed  for  a  measurement,  and  im¬ 
prove  the  accuracy. 
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FIG.  17.  (Color)  Schematic  diagram  of  a  scanning  thermal  microscope 
(SThM).  It  consists  of  a  sharp  temperature- sensing  tip  mounted  on  a  canti¬ 
lever  probe.  The  sample  is  scanned  in  the  lateral  directions  while  the  canti¬ 
lever  deflections  are  monitored  using  a  laser  beam-deflection  technique.  To¬ 
pographical  and  thermal  images  can  be  thermally  obtained.  The  thermal 
transport  at  the  tip-sample  contacts  consists  of  air,  liquid,  and  solid-solid 
conduction  pathways.  A  simple  thermal  resistance  network  model  of  the 
sample  and  probe  combination,  shows  that  when  the  sample  is  at  tempera¬ 
ture  Ts ,  the  tip  temperature  Tt  depends  on  the  values  of  the  thermal  resis¬ 
tances  of  the  tip-sample  contact,  Rts ,  the  tip,  Rt ,  and  the  cantilever  probe, 
Rc. 


B.  Microfabricated  probes 

Access  to  nanofabrication  facilities  offers  new  opportu¬ 
nities  for  probing  directly  the  local  temperature  and  thermal 
transport  mechanisms  in  nanostructures.  Two  significant  ap¬ 
proaches  will  be  commented  on  below:  scanning  thermal  mi¬ 
croscopy  and  the  new  area  of  microinstrumentation  fabrica¬ 
tion. 

A  scanning  thermal  microscope  (SThM)  operates  by 
bringing  a  sharp  temperature- sensing  tip  in  close  proximity 
to  a  sample  solid  surface  (see  Fig.  17).  Localized  heat  trans¬ 
fer  between  the  tip  and  sample  surface  changes  the  tip  tem¬ 
perature.  By  scanning  the  tip  across  the  sample  surface,  a 
spatial  map  of  the  tip- sample  heat  transfer  is  mapped  out. 
When  the  tip  comes  in  local  equilibrium  with  the  sample, 
one  obtains  the  spatial  temperature  distribution  of  the  sample 
surface,  whereas  if  the  temperature  change  is  determined  for 
a  known  heat  flux,  one  could  obtain  the  local  thermal  prop- 
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FIG.  18.  (Color)  Topographical  (top)  and  thermal  images  of  a  multiwall  carbon  nanotube  (MWCN),  10  nm  in  diameter  obtained  by  SThM.  The  bottom  plot 
shows  the  temperature  distribution  along  the  length  of  the  wire  while  the  plot  on  the  side  shows  the  cross-sectional  temperature  profile,  indicating  full-width 
half  maximum  of  the  temperature  peak  to  be  50  nm. 


erties.  The  spatial  resolution  depends  mainly  on  three  factors, 
namely,  tip  sharpness,  tip-sample  heat  transfer  mechanism, 
and  the  thermal  design  of  the  probe.  In  the  recent  past,  the 
spatial  resolution  has  been  improved  to  30-50  nm,  allowing 
one  to  study  thermal  transport  phenomena  at  these  length 
scales.160 

The  idea  of  SThM  was  first  proposed  by  Williams  and 
Wickramasinghe161  when  they  used  a  thermocouple  at  the 
end  of  a  metallic  wire  to  measure  temperature.  Tip-sample 
heat  transfer  was  used  to  image  the  surface  topography  of 
electrically  insulating  materials.  Since  then,  significant 
progress  has  been  made  in  improving  thermocouple-based 
measurements  while  other  techniques  based  on  contact  po¬ 
tential,  electrical  resistance,  and  thermal  expansion  have  also 
been  developed,  which  are  reviewed  in  detail  by 
Majumdar.162  Most  of  the  development  has  focused  on 
cantilever-based  probes  such  that  atomic  force  microscopes 
(AFMs)163  could  be  used  as  a  platform  for  SThM.  When  a 
temperature  sensor  is  mounted  on  the  very  apex  of  the  tip, 
such  probes  can  be  used  to  image  both  the  topography  and 
the  temperature  distribution  of  devices  such  as  single 
transistors26,77  and  vertical  cavity  lasers.164  Such  studies  have 
helped  in  identifying  defects  and  failure  mechanisms  while 
providing  insight  about  electron,  photon  and  phonon  trans¬ 


port  in  these  devices.  In  addition,  SThM  has  also  been  used 
to  measure  thermal  properties  of  materials  and  perform  calo¬ 
rimetry  at  nanometer  scales.165 

There  are  three  important  elements  determining  the  per¬ 
formance  of  SThM,  namely:  (i)  thermal  and  mechanical  de¬ 
sign  of  the  probe;  (ii)  fabrication  of  probe;  (iii)  understand¬ 
ing  of  tip-sample  heat  transfer.  Proper  thermal  design  ensures 
not  only  temperature  measurement  accuracy  but  also  higher 
spatial  resolution.166,167  Fabrication  of  probes  by  MEMS  pro¬ 
cesses  allows  one  to  batch  fabricate  hundreds  of  probes  on  a 
single  wafer  using  optical  lithography,  enabling  multiple  us¬ 
ers  to  utilize  such  probes.160,166  Recent  studies  have  revealed 
the  role  of  gas  conduction,  solid- solid  conduction,  and  con¬ 
duction  through  a  liquid  film  bridging  the  tip  and  the 
sample.164,167  A  fundamental  understanding  of  the  tip-sample 
heat  transfer  allows  accurate  estimation  of  sample  tempera¬ 
ture  distributions  from  the  measured  data.  Recent  progress  in 
these  areas  has  allowed  thermal  imaging  of  nanostructures 
such  as  carbon  nanotubes160  with  30-50  nm  spatial  resolu¬ 
tion  (see  Fig.  18).  Such  a  high  spatial  resolution  allows  one 
to  study  phonon  physics  at  the  length  scale  of  phonon  wave¬ 
lengths  and  mean  free  paths,  which  could  lead  to  greater 
understanding  of  nanoscale  thermal  transport.  In  addition, 
the  spatial  resolution  is  sufficiently  high  to  study  electrother- 
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mal  and  photothermal  effects  in  modern  electronic  and  opto¬ 
electronic  devices. 

Despite  these  recent  advances,  there  are  several  unan¬ 
swered  questions  that  need  to  be  resolved.  It  is  well  known 
that  the  liquid  film  bridging  the  tip  and  the  sample  contrib¬ 
utes  significantly  to  tip-sample  heat  transfer.  Can  the  effect 
of  liquid  conduction  be  controlled  by  tailoring  the  surface 
chemistry?  What  is  the  role  of  near- field  radiation  in  tip- 
sample  heat  transfer,  especially  in  a  vacuum  environment? 
How  do  phonons  and  electrons  propagate  through  a  point 
contact  that  is  on  the  order  of  their  wavelength?  Can  we 
conduct  SThM  experiments  at  low  temperature  and  observe 
quantum  transport  effects?  Because  tip-sample  thermal  resis¬ 
tance  depends  on  surface  topography,  topography-related  ar¬ 
tifacts  in  thermal  images  still  remains  a  problem.  There  are 
two  potential  solutions:  a  thorough  knowledge  of  tip-sample 
heat  transfer  in  order  to  model  the  effect  of  topography;  or  a 
null-point  measurement  such  that  tip-sample  temperature  dif¬ 
ference  can  be  minimized  thereby  eliminating  any  heat  flow. 
The  latter  requires  more  sophisticated  probes  with  multiple 
sensors  and  feedback  control.  Future  progress  will  be  driven 
towards  multiple  sensors  and  actuators  on  integrated  cantile¬ 
ver  probes  such  that  one  could  simultaneously  measure  tem¬ 
perature  and  heat  flux  as  well  as  electrical,  optical,  and  other 
nanoscale  properties  while  imaging  a  nanostructure. 

Finally,  because  new  nanoscale  materials  and  energy 
conversion  devices  are  being  developed,  there  will  also  re¬ 
main  a  need  to  thermally  probe  them  at  nanoscales.  For  ex¬ 
ample,  thermoelectric  properties  of  nanostructured  semicon¬ 
ductors  are  markedly  different  than  those  of  bulk  materials. 
By  simultaneously  measuring  thermal,  electrical,  and  ther¬ 
moelectric  properties  using  a  multifunctional  SThM,  one 
could  characterize  low-dimensional  nanostructures.  An  area 
of  research  that  is  largely  unexplored  is  SThM  under  liquid 
environment,  in  particular,  of  biological  molecules.  Proteins 
and  DNA  undergo  structural  changes  and  phase  transitions  at 
different  temperatures.  Localized  calorimetry  of  molecules 
could  shed  light  on  the  binding  behavior  of  biomolecules. 

One  of  the  major  limitations  of  scanning  probe  micro¬ 
scopes  is  their  inability  to  probe  below  a  solid  surface.  Bal¬ 
listic  electron  emission  microscopy  (BEEM)  allows  one  to 
probe  about  10-50  nm  below  a  conducting  or  a  semiconduct¬ 
ing  surface  by  detecting  electrons  that  travel  ballistically 
through  the  material.168  However,  an  insulating  surface  can¬ 
not  be  probed  by  this  technique.  Because  thermal  transport 
occurs  in  all  materials,  it  may  be  possible  to  use  SThM  to 
probe  materials  below  the  surface  using  either  ballistic  or 
diffusive  phonon  transport.  While  diffusive  phonon  transport 
has  been  explored,165  ballistic  phonon  transport  has  not  yet 
been  studied  using  SThM.  A  combination  of  picosecond  ul¬ 
trasonics  and  SThM  might  perhaps  be  able  to  map  out  ma¬ 
terials  and  their  properties  in  three  dimensions. 

The  SThM  probe  falls  under  a  general  class  of  microin¬ 
struments  that  allows  one  to  thermally  probe  at  a  length  scale 
that  is  an  order  of  magnitude  smaller  than  the  characteristic 
dimension  of  the  instrument.  By  using  traditional  MEMS 
fabrication  techniques,  which  allow  one  to  make  devices  and 
structures  down  to  about  300  nm,  one  could  study  transport 
phenomena  at  length  scales  on  the  order  of  10-50  nm. 
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FIG.  19.  Normalized  differential  reflectivity  for  antimony  at  300  K  (see  Ref. 
170).  Ec  is  the  central  energy  of  the  pulses  which  are  polarized  perpendicu¬ 
lar  to  the  trigonal  axis.  The  data  show  coherent  Alg  phonon  oscillations  of 
frequency  fl0~4.5  THz  and  a  slowly  varying  electronic  background.  The 
fine  structure  near  zero  delay  is  due  to  interference  between  overlapping 
pump  and  probe  beams. 


Hence,  MEMS  provides  a  platform  to  integrate  nanostruc¬ 
tures  and  thereby  forms  a  bridge  between  nano-  and  micros¬ 
cales.  In  the  case  of  the  SThM  probe,  the  nanoscale  tip  is 
integrated  onto  the  rest  of  the  cantilever  and  the  probe.  Fig¬ 
ure  8  shows  another  example  where  two  suspended  micro¬ 
heaters  are  fabricated  on  a  chip  using  standard  MEMS 
techniques.94  A  carbon  nanotube  or  any  ID  nanostructure  can 
be  placed  across  the  two  heater  islands  and  its  thermal  prop¬ 
erties  can  be  measured  using  this  device.  An  advantage  of¬ 
fered  by  MEMS  is  that  hundreds  of  microinstruments  can  be 
built  on  a  single  wafer,  providing  low  cost  measuring  devices 
and  high-throughput  experimentation. 

C.  Coherent  optical  methods 

Coherent  optical  methods,  such  as  stimulated  Raman 
scattering,  have  long  been  applied  to  measurements  of  the 
lifetime  of  optical  phonons.169  As  illustrated  in  Fig.  19,  laser 
excitation  can  lead  to  the  generation  of  intense  vibrations 
showing  a  high  degree  of  temporal  and  spatial  coherence.170 
Transport  measurements  using  optically  generated  coherent 
phonons,  particularly  acoustic  modes,  may  provide  an  alter¬ 
native  to  conventional  thermal  methods  for  the  determination 
of  the  mean  free  path.  Coherent,  monochromatic  phonon 
beams  can  be  generated  using  a  variety  of  techniques  that 
combine  optical  cw  or  pulsed  lasers  and  a  suitable 
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transducer.170-173  Laser-induced  cw  thermomodulation  of  a 
thin  metallic  film  gives  a  tunable  source  of  longitudinally 
polarized  phonons  in  the  1-5  GHz  range.171  Higher  frequen¬ 
cies,  up  to  —0.5  THz,  can  be  attained  using  ultrafast  (subpi¬ 
cosecond)  laser  excitation  of  a  layered  structure.172  In  semi¬ 
conductor  superlattices,  the  mechanism  for  light- sound 
coupling  is  the  acousto-optic  effect,  and  phase  matching  is 
provided  by  the  wave  vectors  involved  in  the  Fourier  decom¬ 
position  of  the  structure.170 

As  discussed  earlier,  it  is  possible  to  generate  coherent 
(but  not  monochromatic)  sound  pulses  as  short  as  100  A 
using  picosecond  ultrasonic  methods.  This  technique  has 
been  recently  applied  to  studies  of  the  phonon  dispersion,173 
and  measurements  of  the  phonon  attenuation  could  be  used 
to  determine  the  mean  free  path.  We  note  that  coherent  pho¬ 
non  sources  can  potentially  be  used  for  microscopy  purposes. 
Here,  there  is  considerable  room  for  improving  the  resolu¬ 
tion.  While  the  wavelengths  of  propagating  modes  in  solids 
are  small  enough  to  resolve  objects  separated  by,  say,  100  A, 
the  length  scale  for  state-of-the-art  ultrasound  and  tomogra¬ 
phy  methods  is  2-3  orders  of  magnitude  larger.174 

VI.  OUTLOOK 

The  many  unresolved  fundamental  issues;  continuing  ad¬ 
vances  in  theory,  computation,  and  experimental  techniques; 
and  important  applications  in  modern  device  materials  en¬ 
sure  that  nanoscale  thermal  transport  will  remain  a  vigorous 
field  of  investigation  for  many  years  to  come.  We  discuss 
specific  examples  of  outstanding  questions  and  recommenda¬ 
tions  for  future  directions  throughout  this  review  but  provide 
a  brief  summary  of  some  selected  examples  here. 

The  limited  range  of  the  data  for  the  thermal  conduc¬ 
tance  of  solid- solid  interfaces  is  not  well  understood — 
theory  predicts  a  much  stronger  dependence  on  acoustic  mis¬ 
match  than  what  is  observed.  Transport  and  coherent  phonon 
experiments  on  individual  and  superlattice  interfaces  with 
systematically  controlled  and  well-characterized  structure  are 
needed.  A  simple,  analytical  model  that  captures  the  basic 
features  of  the  frequency  dependence  in  the  phonon  trans¬ 
mission  coefficient  would  be  of  great  use.  Large-scale  mo¬ 
lecular  dynamics  should  be  used  to  evaluate  the  sensitivity  of 
the  thermal  conductance  to  the  structural  and  chemical  dis¬ 
order  of  the  interface.  A  long  term  goal  is  to  include  elec¬ 
tronic  and  quantum  mechanical  effects  in  atomic-level  simu¬ 
lations  of  thermal  transport. 

Another  important  long  term  challenge  is  the  develop¬ 
ment  of  techniques  that  incorporate  scattering  rates  calcu¬ 
lated  by  atomic-level  simulations  into  higher-level  ap¬ 
proaches,  e.g.,  Monte  Carlo  solutions  of  the  classical  or 
quantum  Boltzmann  transport  equation.  Improvements  in 
scanning  thermal  microscopy  and  other  nanoscale  metrology 
methods  are  needed  to  make  quantitative  comparisons  be¬ 
tween  experimental  observations  of  nanoscale  transport  and 
the  predictions  of  theory. 
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